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INTRODUCTION: 


A  new  concept  of  flat-panel  imager  (FPI)  with  avalanche  gain  and  high  resolution  (with  50  micron 
pixel  size)  is  being  investigated  for  improving  the  imaging  performance  of  digital  mammography  at 
low  dose  and  high  spatial  frequencies,  which  are  critical  for  the  detection  of  subtle  breast  abnonnalities 
and  the  development  of  digital  tomosynthesis.  The  detector  employs  an  avalanche  photoconductor  - 
amorphous  selenium  (a-Se),  which  is  called  HARP,  to  detect  and  amplify  the  optical  signal  generated 
by  an  x-ray  scintillator  with  needle  structure  -  cesium  iodide  (Csl),  and  form  a  charge  image  that  is 
read  out  by  a  high-resolution  field  emitter  array  (FEA).  We  call  the  proposed  detector  SAPHIRE 
(Scintillator-Avalanche  Photoconductor  with  High  Resolution  Emitter  readout).  Our  objective  is  to 
investigate  the  feasibility  of  SAPHIRE  and  demonstrate  its  potential  advantages  over  existing  detector 
technology. 

BODY: 

The  following  tasks  were  proposed  in  the  approved  Statement  of  Work.  The  results  of  our  investigation 
in  the  first  year  are  summarized  as  follows: 

TASK  1:  Detennine  the  optimal  CsI-HARP  combination  for  mammography,  Month  1-24 

a.  Develop  a  linear  system  model  for  SAPHIRE,  so  that  the  imaging  performance  of  different 
CsI-HARP  combination  can  be  predicted.  Month  1-6 

In  order  to  detennine  the  optimal  Csl  layer  for  SAPHIRE,  we  investigated  the  imaging  performance  of 
four  different  samples  of  Csl  layers  with  thickness  of  150  pm,  which  is  appropriate  for  mammography. 
The  Csl  samples  differ  in  theirs  manufacturing  methods.  Two  were  optimized  for  high  resolution  (HR) 
and  the  other  two  were  optimized  for  high  light  (HL).  For  each  type  (HR  or  HL)  of  Csl,  two  different 
substrates  were  used.  The  samples  with  fiber  optic  face  plate  were  called  FOS,  and  the  samples  with 
amorphous  carbon  substrates  were  called  ACS.  We  found  that  ACS  HR  plate  has  higher  light  output 
than  the  FOS  HR  plate,  however  it  has  slightly  lower  Swank  factor  (i.e.  higher  added  noise)  due  to  the 
higher  variance  in  light  output  for  x-rays  absorbed  at  different  depth  of  the  Csl  layer.  This  is  directly 
related  to  the  method  of  manufacturing  of  the  Csl  samples.  This  investigation  has  generated  a  lot  of 
interests  from  the  manufacturer,  Hamamatsu  Photonics,  and  we  will  be  collaborating  on  the 
optimization  of  Csl  manufacturing  methods  for  different  x-ray  imaging  applications.  This  work  is 
described  in  detail  in  a  Med.  Phys.  paper  published  earlier  this  year  (manuscript  1). 

b.  Make  a  test  system  for  experimental  measurement  of  the  imaging  performance  of  CsI- 
HARP  combination  using  a  HARP  tube  with  fiber  optic  faceplate.  Month  1-12 

We  made  an  improved  HARP  tube  with  FOP  faceplate.  It  has  much  fewer  defects  than  the  first  one 
made  last  year,  the  thickness  of  the  HARP  layer  was  increased  to  7  pm  (from  the  4  pm  last  year). 
Shown  in  Figure  1  (left)  is  a  picture  of  the  FOP  HARP  tube  installed  in  an  NTSC  standard  camera  with 
separate  high  voltage  control  for  the  target. 
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Figure  1:  Left:  Improved  prototype  FOP  FIARP  tube  installed  in  a  NTSC  camera.  The  FOP  face  plate 
is  in  direct  contact  with  a  Csl  layer  at  the  front  surface.  Right:  a  flat  field  x-ray  image  made  with  the 
CsI-FOP  HARP  combination,  point  defects  are  visible. 


Figure  2:  X-ray  images  obtained  with  the  prototype  FOP  HARP  directly  coupled  to  a  Csl  layer  (150 
pm  HL  type).  Left:  image  of  a  dog  heart  and  a  180  pm  thick  wire;  Right:  x-ray  image  of  a  bar  pattern. 
The  resolution  of  the  bars  from  right  to  left  are:  6.3,  8,  10,  12.5,  16  and  20  lp/mm,  respectively.  It 
shows  that  the  bars  with  16  lp/mm  are  visible. 

c.  Measure  the  x-ray  imaging  performance  of  Csl -HARP  using  Csl  of  different  types  and 
thickness,  and  determine  the  optimal  CsI-HARP  combination.  Month  13-24 

A  150  pm  thick  Csl  layer  (HL  type  optimized  for  maximum  light  output)  was  coupled  directly  to  the 
FOP  faceplate  for  x-ray  detection.  Shown  in  Figure  1  (right)  is  a  flat-field  x-ray  image  made  with  the 
Csl  coupled  FOP  HARP.  The  gain  of  the  HARP  tube  is  approximately  60.  Only  several  point  defects 
are  visible  on  the  right  hand  side  of  the  active  imaging  area,  which  is  a  significant  improvement  over 
the  large  area  of  defects  of  the  previous  prototype.  This  was  made  possible  by  mechanical  polishing  of 
the  FOP  prior  to  soft  (optical)  polishing  before  HARP  deposition.  This  prototype  FOP  HARP  tube 
would  allow  us  to  conduct  systematic  characterization  of  the  imaging  properties  of  the  CsI-HARP 
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combination.  Figure  2  (left)  shows  is  an  x-ray  image  of  a  dog  heart  with  a  180  gm  wire,  and  Figure  2 
(right)  is  the  image  of  a  resolution  bar  pattern.  It  shows  that  while  the  FOP  coupling  between  Csl  and 
HARP  is  able  to  preserve  the  resolution,  there  is  some  image  blur  mostly  caused  by  the  inherent  image 
resolution  of  the  Csl  layer,  which  is  optimized  for  light  output.  These  preliminary  results  are 
encouraging.  We  are  going  to  start  a  full  characterization  of  the  MTF,  NPS  and  DQE  of  different  Csl 
layers  in  combination  with  HARP,  and  determine  the  advantage  of  avalanche  gain,  as  predicted  in  our 
theoretical  investigation. 

TASK  2:  Investigate  the  factors  affecting  the  spatial  resolution  of  the  FEA  readout  method,  Month  1- 
36 


a.  Develop  an  electric  field  simulation  for  the  electron  optics  of  SAPHIRE  to  calculate  the 
electron  beam  trajectory  for  FEA  readout  method.  Month  1-12 

b.  Determine  the  optimal  mesh  electrode  configuration  that  results  in  the  highest  spatial 
resolution.  Month  13-18 

c.  Refine  optimal  mesh  electrode  configuration  in  connection  with  the  noise  model.  Month 
19-21 

d.  Calculate  the  spatial  resolution  due  to  FEA  readout  for  the  prototype  detector  to  verily 
experimental  measurement.  Month  25-36 

Our  previous  studies  showed  that  by  placing  the  HARP  target  and  the  FEA  close  to  each  other  (0.5 
mm)  and  applying  a  high  mesh  voltage  that  is  >400  V,  the  lateral  spreading  of  the  electron  beam 
emitted  from  the  FEA  could  be  limited  to  <  50  pm.  However  it  is  very  challenging  in  practice  to 
maintain  a  vacuum  spacing  of  0.5  over  a  large  area,  practical  distance  is  usually  -  1  mm.  Furthennore, 
the  electron  beam  spread  is  a  function  of  signal  charge,  which  makes  the  images  non-uniform. 

Therefore  we  conducted  further  investigation  of  the  methods  for  improving  spatial  resolution  of 
SAPHIRE.  Our  simulation  showed  that  by  adding  an  integrated  focusing  electrode  on  the  FEA,  the 
electron  beam  spread  can  be  limited  to  within  2  pm,  which  makes  the  size  of  the  electron  beam  defined 
by  the  pixel  size  of  the  FEA  instead  of  the  lateral  spread.  Shown  in  Figure  3  are  the  results  of  our 
simulation,  which  is  the  trajectory  of  electrons  emitted  from  the  FEA  (y=0)  before  reaching  the  HARP 
target. (y=  1000  pm).  The  distance  between  HARP  and  FEA  was  1  mm,  with  a  mesh  voltage  of  350  V. 

It  shows  that  without  focusing  electrode,  the  lateral  spread  of  the  electron  beam  emitted  from  each 
FEA  tip  is  ~  60  pm  (per  side),  which  is  not  appropriate  for  pixel  sizes  <  100  pm.  With  the  focusing 
electrode,  which  is  an  integrated  electrode  held  at  negative  potential  above  the  gate  electrode  of  the 
FEA,  the  lateral  spread  of  the  electron  beam  is  significantly  reduced. 

A  manuscript  is  currently  being  prepared  for  submission  entitled:  “SAPHIRE  (Scintillator-Avalanche 
Photoconductor  with  High  Resolution  Emitter  readout)  for  X-ray  Imaging:  Spatial  Resolution  and 
Lag”.  It  includes  our  theoretical  investigation  of  the  factors  affecting  the  spatial  resolution  and  image 
lag  of  SAPHIRE. 
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Lateral  Spread  (|am) 

Figure  3:  Simulated  electron  trajectory  for  FEA-HARP  distance  of  1  mm.  The  mesh  is  placed 

in  the  center,  i.e.  0.5  mm  from  the  FEA  and  HARP,  and  held  at  350  V  (not  shown  in  the  graph). 

The  lateral  spread  of  electrons  exceeds  60  pm  on  each  side.  With  the  focusing  electrode  held  at 

a  negative  potential  VL,  the  beam  spread  depends  can  be  reduced  to  <  2pm. 

TASK  3:  Investigate  the  factors  affecting  the  electronic  noise  and  readout  speed,  Month  1-21 

a.  Investigate  the  factors  affecting  the  readout  speed  of  SAPHIRE  and  detennine  the  readout 
time  as  a  function  of  the  number  of  ITO  stripes  and  characteristics  of  charge  amplifier 
ASIC.  Month  7-12 

b.  Investigate  the  factors  affecting  the  electronic  noise  and  develop  a  noise  model  based  on  the 
properties  of  charge  amplifiers,  the  number  of  divided  ITO  stripes,  the  readout  speed  and 
the  electron  optics  design.  Month  13-18 

c.  Optimize  the  number  of  stripes  and  the  mesh  electrode  configuration  so  that  the  best  noise 
performance  and  adequate  readout  speed  can  be  achieved.  Month  19-21 

We  detennined  that  in  order  to  keep  the  electronic  noise  below  2000  electrons,  the  ITO  signal 
electrode  for  HARP  should  be  divided  into  strips.  For  a  full  field  digital  mammography  detector  of  24 
cm  x  30  cm,  dividing  into  32  strips  is  sufficient  to  keep  the  amplifier  input  load  capacitance  <  50  pF, 
which  leads  to  an  amplifier  noise  of  <  2000  electrons.  Division  of  the  signal  electrode  also  permits 
parallel  readout  of  signal  by  firing  32  FEA  pixels  at  the  same  time.  Although  readout  speed  in 
mammography  is  not  demanding  at  30  second  per  image,  the  new  development  of  digital 
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tomosynthesis  requires  a  frame  rate  of  3  frames  /  second.  The  readout  lag  of  SAPHIRE  will  be 
dominated  by  the  beam  discharge  lag  of  the  electrons,  which  is  a  function  of  the  energy  spread  of  the 
electron  beam  and  the  thickness  of  HARP.  Figure  4  shows  the  first  frame  beam  discharge  lag  as  a 
function  of  frame  time  for  two  different  target  thickness  values.  We  assumed  that  the  ITO  electrode  of 
HARP  is  divided  into  32  strips,  and  the  pixel  size  is  50  pm.  This  means  that  the  detector  has  6000  x 
4800  pixels,  and  each  strip  serves  as  the  signal  electrode  for  6000  x  150  pixels.  Figure  4  shows  that  for 
a  HARP  thickness  of  8  pm  and  a  frame  rate  of  3  frames  per  second,  the  first  frame  lag  (at  300  ms)  is 
less  than  2  %.  For  a  thicker  target  with  smaller  capacitance,  the  lag  is  negligible. 

TASK  4:  Demonstrate  the  practical  feasibility  of  SAPHIRE  by  building  a  small  area  prototype 
detector,  Month  22-36 

a.  Construct  the  small  area  prototype  detector  according  to  the  optimal  mesh  electrode 
configuration  determined  from  TASK  3.  Month  22-24 

b.  Evaluate  the  optical  imaging  perfonnance  of  the  prototype  detector  without  Csl  at  the  input. 
Month  25-27 

c.  Evaluate  the  x-ray  imaging  performance  of  the  prototype  detector  with  Csl  coupled  to  the 
input.  Month  28-36 


KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Developed  a  cascaded  linear  system  model  for  SAPHIRE. 

•  Investigated  the  imaging  performance  of  Csl  layers  made  with  different  manufacturing 
techniques,  and  determined  that  the  Csl  optimized  for  high  resolution  (HR)  can  provide  better 
DQE  at  higher  spatial  frequencies  required  for  mammography.  It  was  also  discovered  that  HR 
Csl  deposited  on  FOP  substrates  provides  the  highest  Swank  factor,  whereas  HR  Csl  on 
amorphous  carbon  substrate  provides  higher  light  output  at  the  cost  of  a  lower  Swank  factor. 

•  Investigated  the  factors  affecting  the  spatial  resolution  and  image  lag.  Determined  that  dividing 
the  ITO  signal  electrode  of  a  full  field  digital  mammography  detector  into  32  strips  is  sufficient 
for  a  readout  speed  of  3  frames  per  second,  and  keeping  the  electronic  noise  to  <  2000 
electrons. 

•  Constructed  a  prototype  HARP  tube  with  FOP  substrate,  which  can  be  coupled  directly  to 
different  types  of  Csl  layers. 


REPORTABLE  OUTCOME: 

Manuscripts: 

1.  Wei  Zhao,  Dan  Li,  Alla  Reznik,  B.  J.  M.  Lui,  D.  C.  Hunt,  J.  A.  Rowlands,  Yuji  Ohkawa  and 
Kenkichi  Tanioka,  “Indirect  flat-panel  detector  with  avalanche  gain:  Fundamental  feasibility 
investigation  for  SHARP-AMFPI  (Scintillator  HARP  Active  Matrix  Flat  Panel  Imager)”,  Med. 
Phys.  2954-2966  (2005) 
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2.  A.  R.  Lubinsky,  Wei  Zhao,  Goran  Ristic  and  J.  A.  Rowlands,  “Screen  optics  effects  on  detective 
quantum  efficiency  in  digital  radiography:  Zero-frequency  effects”,  Med.  Phys.  33,  1499-1509 
(2006) 

3.  A.  Reznik,  B.  J.  M.  Lui,  Y.  Ohkawa,  T.  Matsubara,  K.  Miyakawa,  M.  Kubota,  K.  Tanioka,  T. 
Kawai,  W.  Zhao  and  J.  A.  Rowlands,  “The  quantum  effience  of  photo-charge  generation  in  a-Se 
avalanche  photodetectors”,  Nuclear  Instrument  and  Methods  A  (accepted) 

Conference  presentations  and  published  abstract: 

1  Dan  Li,  Wei  Zhao,  Geordi  Pang  and  J.  A.  Rowlands,  “A  new  concept  of  indirect  flat-panel 
detector:  SAPHIRE  (Scintillator  Avalanche  Photoconductor  with  High  Resolution  Emitter 
readout)”,  AAPM  Annual  Conference,  August,  2005,  Seattle,  WA 


CONCLUSIONS 

During  the  second  year  of  our  research  supported  by  Army  BCRP,  we  have  focused  our  investigation 
on  understanding  the  factors  affecting  the  spatial  resolution,  image  lag  and  electronic  noise  of 
SAPHIRE.  We  found  that  adding  a  focusing  electrode  at  each  FEA  tip  can  reduce  the  lateral  spreading 
of  the  electron  beam  from  60  pm  to  2  pm.  The  focusing  electrode  allows  the  distance  between  HARP 
and  the  FEA  to  be  kept  at  a  practical  distance,  i.e.  1  mm  or  more,  and  the  potential  of  the  mesh 
electrode  reasonably  low,  e.g.  350  V,  while  maintaining  a  high  resolution  that  is  defined  by  the  pixel 
size  of  the  FEA  array.  By  dividing  the  ITO  signal  electrode  to  a  moderate  number  of  strips,  e.g.  32,  the 
electronic  noise  can  be  kept  below  2000  electrons,  and  the  readout  rate  of  a  full  field  digital 
mammography  detector  with  50  pm  pixel  size  could  be  3  frames  per  second  or  higher. 

REFERENCES 
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Indirect  flat-panel  detector  with  avalanche  gain:  Fundamental  feasibility 
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Alla  Reznik,  B.  J.  M.  Lui,  D.  C.  Hunt,  and  J.  A.  Rowlands 

Imaging  Research,  Sunnybrook  and  Women ’s  Health  Sciences  Center,  2075  Bayview  Avenue,  Toronto, 

Ontario,  M4N  3M5  Canada 

Yuji  Ohkawa  and  Kenkichi  Tanioka 

Advanced  Imaging  Devices  Research  Division,  Science  &  Technical  Research  Laboratories, 

Japan  Broadcasting  Corporation,  1-10-11  Kinuta,  Setagaya-ku,  Tokyo  157-8510  Japan 

(Received  11  April  2005;  revised  30  June  2005;  accepted  for  publication  30  June  2005; 
published  30  August  2005) 

An  indirect  flat-panel  imager  (FPI)  with  avalanche  gain  is  being  investigated  for  low-dose  x-ray 
imaging.  It  is  made  by  optically  coupling  a  structured  x-ray  scintillator  CsI(Tl)  to  an  amorphous 
selenium  (a-Se)  avalanche  photoconductor  called  HARP  (high-gain  avalanche  rushing  photocon¬ 
ductor).  The  final  electronic  image  is  read  out  using  an  active  matrix  array  of  thin  film  transistors 
(TFT).  We  call  the  proposed  detector  SHARP-AMFPI  (scintillator  HARP  active  matrix  flat  panel 
imager).  The  advantage  of  the  SHARP-AMFPI  is  its  programmable  gain,  which  can  be  turned  on 
during  low  dose  fluoroscopy  to  overcome  electronic  noise,  and  turned  off  during  high  dose  radi¬ 
ography  to  avoid  pixel  saturation.  The  purpose  of  this  paper  is  to  investigate  the  important  design 
considerations  for  SHARP-AMFPI  such  as  avalanche  gain,  which  depends  on  both  the  thickness  c/Se 
and  the  applied  electric  field  ESe  of  the  HARP  layer.  To  determine  the  optimal  design  parameter  and 
operational  conditions  for  HARP,  we  measured  the  ESe  dependence  of  both  avalanche  gain  and 
optical  quantum  efficiency  of  an  8  p,m  HARP  layer.  The  results  were  used  in  a  physical  model  of 
HARP  as  well  as  a  linear  cascaded  model  of  the  FPI  to  determine  the  following  x-ray  imaging 
properties  in  both  the  avalanche  and  nonavalanche  modes  as  a  function  of  ESe:  (1)  total  gain  (which 
is  the  product  of  avalanche  gain  and  optical  quantum  efficiency);  (2)  linearity;  (3)  dynamic  range; 

(4)  gain  nonuniformity  resulting  from  thickness  nonuniformity;  and  (5)  effects  of  direct  x-ray 
interaction  in  HARP.  Our  results  showed  that  a  HARP  layer  thickness  of  8  /rm  can  provide  ad¬ 
equate  avalanche  gain  and  sufficient  dynamic  range  for  x-ray  imaging  applications  to  permit  quan¬ 
tum  limited  operation  over  the  range  of  exposures  needed  for  radiography  and 
fluoroscopy.  ©  2005  American  Association  of  Physicists  in  Medicine. 

[DOI:  10.1118/1.2008428] 
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I.  INTRODUCTION 

An  active  matrix  flat-panel  imager  (AMFPI)  based  on  thin 
film  transistor  (TFT)  arrays  is  the  most  promising  technology 
for  digital  x-ray  imaging  due  to  its  compact  size,  rapid  read¬ 
out,  and  superior  performance  compared  to  screen  films.1 
AMFPI  are  categorized  as  either  direct  or  indirect  based  on 
the  materials  used  for  x-ray  detection  (x-ray  photoconductors 
or  scintillators,  respectively).  Direct  AMFPI  has  the  advan¬ 
tages  of  higher  image  resolution  and  a  simpler  TFT  array 
structure  that  can  be  manufactured  in  a  standard  facility  for 
active  matrix  liquid  crystal  displays  (AMLCD).  Despite  their 
differences,  both  types  of  AMFPI  offer  better  image  quality 
than  screen  films  and  computed  radiography  (CR).2'3 

However  before  AMFPI  can  equal  the  performance  of  the 
x-ray  image  intensifier  (XRII)  in  fluoroscopy,  two  major  dif¬ 
ficulties  must  be  overcome:  (i)  the  ability  to  produce  good 


image  quality  at  a  very  low  dose,  such  as  in  the  dark  part  of 
a  fluoroscopy  image  (—0.1  //R  per  frame)4’3  and  (ii)  the  abil¬ 
ity  to  produce  images  at  a  high  frame  rate  without  artifacts, 
especially  when  the  radiation  exposure  is  switched  from  ra¬ 
diographic  to  fluoroscopic.6’7  Even  in  well-established  radio- 
graphic  applications,  AMFPI  are  facing  challenges  because 
their  rapid  readout  and  freedom  from  geometric  distortion 
have  generated  a  renewed  interest  in  advanced  clinical  appli- 
cations.  One  example  is  digital  tomosynthesis,  in  which  a 
rapid  sequence  of  images  is  taken  from  different  angles,  and 
tomographic  slices  at  different  depths  of  the  patient  anatomy 
are  reconstructed.  The  application  of  AMFPI  in  cone -beam 
computed  tomography9'10  (CT)  is  also  under  intensive  inves¬ 
tigation.  Such  advanced  applications  of  AMFPI  pose  chal¬ 
lenges  for  the  development  of  the  next  generation  of  detec¬ 
tors,  primarily  regarding  the  ability  of  AMFPI  to  generate 
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high  quality  images  that  are  quantum  noise  limited  and  free 
from  artifacts  at  low  x-ray  exposures  and  high  frame  rates.  It 
has  been  shown  that  the  detective  quantum  efficiency  (DQE) 
of  AMFPI  is  inferior  to  XRII  at  the  lower  end  of  the  dose 
range  (0.1-1  //R  per  frame)  encountered  in  fluoroscopy.11  In 
mammography,  the  detector  is  not  x-ray  quantum  noise  lim¬ 
ited  behind  dense  breast,  where  the  average  exposure  is 
~1  mR,  and  the  DQE  of  the  detector  is  severely  compro- 

12  13 

mised  by  the  electronic  noise.  ’ 

The  strategies  for  improving  the  low  dose  performance  of 
AMFPI  can  be  divided  into  two  categories:  one  is  to  increase 
the  x-ray  to  image  charge  conversion  gain  so  that  the  x-ray 
quantum  noise  can  overcome  the  electronic  noise  and  the 
other  is  to  decrease  the  electronic  noise.  These  strategies  are 
common  to  both  types  of  AMFPI  since  they  have  approxi¬ 
mately  the  same  conversion  gain  and  pixel  electronic  noise.1 
Several  investigations  are  taking  place  to  increase  the  gain 
for  direct  AMFPI  by  developing  other  x-ray  photoconduc¬ 
tors.  Promising  results  have  been  shown  for  lead  iodide 
(Pbl2)  and  mercuric  iodide  (Hgl2).14  The  conversion  gain  of 
these  two  materials  is  5-7  times  higher  than  that  of  current 
systems.  The  challenges  for  this  technology  are  to  develop 
deposition  techniques  for  large  area  detectors  with  good  im¬ 
aging  properties,  and  to  ensure  that  both  electrons  and  holes 
have  adequate  range  to  reduce  the  probability  of  charge 
trapping.13  To  reduce  electronic  noise,  several  investigators 
have  proposed  advanced  pixel  designs,  which  incorporate 
pixel  amplification  by  adding  at  least  two  more  TFTs  at  each 
pixel. 16-18  The  electronic  noise  can  thus  be  reduced  to  ~700 
electrons  per  pixel,  which  is  a  factor  of  2-4  lower  than  ex¬ 
isting  AMFPI.  A  challenge  for  this  approach  is  the  practical 
implementation  of  the  complex  pixel  design  over  a  large  area 
with  consistent  and  uniform  imaging  performance  because 
each  pixel  operates  as  an  analog  amplifier  circuit  as  opposed 
to  the  simple  switching  device  (digital)  in  existing  AMFPI. 
Due  to  these  practical  challenges  none  of  the  above  ap¬ 
proaches  for  increased  gain  and  reduced  electronic  noise  has 
been  implemented  in  commercial  large  area  AMFPI. 

We  are  proposing  a  new  concept  of  AMFPI  by  adding 
avalanche  gain  to  an  indirect  AMFPI,  as  shown  in  Fig.  1.  The 
optical  photons  emitted  from  the  Csl  layer  upon  absorption 
of  x  rays  are  detected  and  amplified  by  an  avalanche  a-Se 
photoconductor  called  HARP  (high-gain  avalanche  rushing 
photoconductor).  Light  photons  absorbed  by  the  HARP  layer 
generate  electron-hole  pairs  near  the  top  interface 
(<0.1  /iin).  Under  a  sufficiently  high  ESe  (produced  by  ap¬ 
plying  positive  voltage  on  the  top  bias  electrode  of  a-Se), 
holes  moving  toward  the  bottom  surface  will  undergo  ava¬ 
lanche  multiplication  and  create  more  holes  and  electrons. 
The  avalanche  gain  gav  ranges  from  unity  up  to  ~  1000  times 
depending  on  ESe  and  the  thickness  of  the  a-Se  layer,  dSe.19 
The  electrons  move  to  the  top  bias  electrode  while  holes  are 
collected  by  the  pixel  electrodes  and  form  a  charge  image, 
which  is  readout  with  a  two-dimensional  array  of  TFTs,  simi¬ 
lar  to  that  in  existing  AMFPI.  The  proposed  detector  is 
named  SHARP-AMFPI  (scintillator-HARP  active  matrix 
flat-panel  imager).  Rougeot  and  Possin  have  suggested  a 
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Fig.  1 .  (a)  Diagram  showing  the  concept  of  SHARP-AMFPI,  where  an  ava¬ 
lanche  photoconductor  is  used  to  detect  light  photons  generated  by  Csl.  (b) 
Side  view  showing  the  multilayer  structure  of  the  detector  (thickness  not 
drawn  to  scale). 
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similar  concept,  but  it  has  not  been  implemented.  The  only 

other  indirect  AMFPI  approach  involving  an  a-Se  photocon- 

2 1 

ductor  as  an  optical  sensor  did  not  have  avalanche  gain. 

Compared  to  existing  indirect  FPI,  SHARP-AMFPI  has 
the  following  advantages:  (1)  x-ray  quantum  noise  limited  at 
a  much  lower  x-ray  exposure  level;  (2)  better  temporal  per¬ 
formance,  i.e.,  lower  lag  and  ghosting  because  of  the  lower 
trap  density  in  a-Se  (compared  to  a-Si)  and  the  much  higher 
electric  field  ESe;  (3)  compatibility  with  standard  TFT  array 
manufacturing  facility  for  AMLCD;  (4)  ability  to  use  high- 
resolution  (HR)  type  Csl,  which  cannot  be  used  in  existing 
AMFPI  due  to  its  lower  light  output.  Compared  to  the  high 
conversion  gain  x-ray  photoconductors  being  developed  for 
direct  FPI,  the  proposed  detector  has  the  following  advan¬ 
tages:  (1)  Programmable  gain  by  changing  the  electric  field 
ESe  as  opposed  to  the  fixed  high  conversion  gain  of  Hgl2  or 
Pbl2.  One  of  the  practical  problems  of  having  a  fixed  high 
gain  is  that  the  large  signal  charge  (at  high  exposure)  re¬ 
quires  a  large  pixel  storage  capacitor  (>15  pF),  which  is 
impractical  to  make  especially  for  small  pixel  sizes.  This 
means  that  a  detector  that  is  optimized  for  low  dose  (fluoro¬ 
scopic)  x-ray  imaging  applications  will  not  work  properly  at 
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high  radiographic  dose  levels.  With  SHARP-AMFPI,  by  ap¬ 
propriately  decreasing  the  ESe  during  high  exposure  applica¬ 
tions,  the  signal  charge  is  comparable  to  that  at  low  dose 
with  avalanche  gain.  This  simplifies  the  TFT  array  pixel  de¬ 
sign  and  permits  a  detector  with  wide  dynamic  range.  (2) 
Better  temporal  performance  because  both  charge  carriers  in 
a-Se  have  adequate  range  and  the  probability  of  charge  trap¬ 
ping  is  low.  This  is  in  contrast  to  the  very  short  range  of  one 
type  of  carriers  in  Pbl2  (electrons)  and  Hgl2  (holes)  which 
results  in  significant  charge  trapping  and  ghosting.  The  ad¬ 
vantages  described  above  make  the  proposed  detector  suit¬ 
able  for  all  x-ray  imaging  modalities  (i.e.,  radiography  and 
fluoroscopy).  The  purpose  of  this  paper  is  to  investigate  the 
important  design  parameters  for  HARP  and  to  establish  the 
operational  conditions  that  ensure  adequate  dynamic  range 
for  both  fluoroscopy  and  radiography. 

Since  avalanche  multiplication  in  a- Se  was  discovered 
over  a  decade  ago  the  HARP  video  tube,  which  consists  of  a 
HARP  layer  readout  with  a  scanning  electron  beam,  has  been 
commercialized  for  high  sensitivity  and  high  definition  tele¬ 
vision  applications.  Stable  and  uniform  avalanche  multipli- 

22 

cation  has  been  demonstrated  with  very  little  added  noise." 
No  correction  for  fixed  pattern  noise  is  necessary.  The  elec¬ 
tric  field  ESe  under  which  avalanche  multiplication  occurs  is 
>70  V/um,  which  is  an  order  of  magnitude  higher  than  that 
typically  used  in  direct  AMFPI  incorporating  thick  a-Se  lay¬ 
ers.  To  sustain  Es e  of  this  high  magnitude  without  a  signifi¬ 
cant  increase  in  dark  current  due  to  charge  injection  from  the 
bias  electrodes,  blocking  layers  are  required  for  both  elec¬ 
trode  interfaces.  As  shown  in  Fig.  1(b),  the  HARP  sensor  is  a 
multilayer  structure.19  It  contains  a  transparent  indium  tin 
oxide  (ITO)  electrode  which  is  biased  positively  during  im¬ 
aging,  followed  by  a  thin  layer  (—20  nm)  of  cerium  oxide 
(Ce02)  and  a  thin  layer  of  lithium  fluoride  (LiF)  doped  a-Se 
to  block  injection  of  holes  from  the  ITO.  The  intrinsic  a-Se 
photoconductor  is  from  0.5  to  35  /rm  thick  depending  on  the 
desired  avalanche  gain.  The  bottom  interface  of  the  HARP 
structure  is  a  submicron  thin  layer  of  antimony  trisulfide 
(Sb2S3),  which  is  used  to  block  electron  injection  from  the 
negative  bias  electrode  (i.e.,  pixel  electrodes  of  the  TFT  ar¬ 
ray).  Since  the  mobility  of  holes  in  a-Se  is  more  than  30 
times  higher  than  that  of  electrons,  only  holes  create  a  sig¬ 
nificant  number  of  carriers  by  impact  ionization.  Hence  it  is 
essential  for  light  photons  to  enter  from  the  positive  bias 
electrode. 

Our  previous  calculation  using  a  preliminary  linear  sys¬ 
tem  model  of  indirect  FPI  with  avalanche  gain  predicted  that 
an  avalanche  gain  of  gav=20  is  needed  to  produce  x-ray 
quantum  noise  limited  performance  for  fluoroscopy  and  low 
dose  mammography  (tomosynthesis).23  The  avalanche  gain 
of  HARP  increases  as  a  function  of  both  the  electric  field  ESe 
and  the  thickness  t/Se.  To  minimize  dark  current  ESe  is  usu¬ 
ally  1 20  V/  j-L m .  A  value  of  gav=20  requires  a  HARP  layer 
thickness  of  4  fim  or  above.  For  the  present  paper  we  in¬ 
vestigated  the  effects  of  the  optical  imaging  properties  of 
HARP  on  the  x-ray  imaging  performance  of  SHARP- 
AMFPI,  so  that  appropriate  design  parameters  and  opera- 


Fig.  2.  Experimental  apparatus  for  measuring  the  photosensitivity  of  HARP. 

tional  conditions  could  be  determined.  We  measured  the  ESe 
dependence  of  both  avalanche  gain  and  optical  quantum  ef¬ 
ficiency  of  an  8  fim  HARP  layer.  The  results  were  used  in  a 
physical  model  of  HARP  as  well  as  a  linear  cascaded  model 
of  the  FPI  to  determine  the  following  x-ray  imaging  proper¬ 
ties  in  both  the  avalanche  and  nonavalanche  modes  as  a 
function  of  dSe  and  ESe:  (1)  total  gain  (which  is  the  product 
of  avalanche  gain  and  optical  quantum  efficiency);  (2)  lin¬ 
earity;  (3)  dynamic  range;  (4)  gain  nonuniformity  resulting 
from  thickness  nonuniformity;  and  (5)  the  effects  of  direct 
x-ray  interaction  in  HARP. 

II.  PHOTOSENSITIVITY  OF  HARP 

The  x-ray  imaging  performance  of  indirect  FPI  with  ava¬ 
lanche  gain  relies  on  the  photosensitivity  of  HARP.  Under¬ 
standing  the  dependence  of  photosensitivity  on  ESe  and  the 
layer  thickness  <7Se  is  crucial  for  guiding  the  choice  of  HARP 
design  parameters  and  operational  conditions  for  both  the 
avalanche  and  nonavalanche  modes.  In  this  section,  we  will 
first  present  the  experimental  method  and  results  of  the 
HARP  photosensitivity  measurements,  and  then  describe  the 
theory  of  ESe  dependence  of  photosensitivity,  which  is  pro¬ 
portional  to  both  the  optical  quantum  efficiency  (QE)  rj  and 
the  avalanche  gain  gav.  The  theory  and  experimental  data  will 
then  be  used  to  determine  77  and  gav  as  a  function  of  ESe  in 
the  context  of  x-ray  imaging. 

A.  Experimental  measurement  of  photosensitivity 
1.  Method 

The  experimental  apparatus  for  photosensitivity  measure¬ 
ment  is  shown  in  Fig.  2.  A  grating  monochromator  was  used 
to  generate  optical  photons  of  a  single  wavelength  The 
values  of  A  used  in  our  experiments  were  in  the  range  of 
380-600  nm.  The  intensity  of  the  output  beam  was  attenu¬ 
ated  by  neutral-density  (ND)  filters.  A  beamsplitter  was  used 
to  direct  a  fraction  (90%)  of  the  attenuated  beam  to  a  silicon 
(Si)  photodiode  (PD),  which  was  used  to  monitor  the  beam 
intensity  during  the  experiment.  The  optical  photons  passing 
directly  through  the  beam  splitter  (10%)  were  detected  by  the 
HARP  video  tube.  Before  photosensitivity  measurements  the 
input  light  intensity  to  the  HARP  was  calibrated  such  that  a 
constant  input  power  was  used  for  each  X.  This  was  achieved 
by  adjustment  of  the  ND  filters  based  on  the  light  intensity 
monitored  by  the  Si  PD. 

For  each  A.  and  ESe,  the  signal  current  S  from  the  HARP 
tube  was  recorded.  The  photosensitivity  was  quantified  as  the 
number  of  signal  charge  carriers  generated  per  absorbed  light 
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Fig.  3.  Measured  effective  quantum  efficiency  if  as  a  function  of  ESe  for  an 
8  /nm  thick  HARP  layer. 


photon.  This  is  referred  to  as  the  effective  quantum  effi¬ 
ciency  rf ,  and  it  is  the  product  of  77  and  gav.  The  values  of  r; 
were  experimentally  determined  from  the  measurement  of 
HARP  tube  current  S  using 


Sle 

ITJhv 


(1) 


where  e  is  the  electron  charge,  1  is  the  HARP  input  light 
intensity  power  given  in  Watts  (W),  hv  is  the  energy  of  each 
incident  photon,  and  Ts  is  the  fraction  of  the  photons  incident 
on  the  HARP  tube  reaching  the  a-Se  layer.  The  correction 
factor  Ts  represents  the  light  attenuation  by  the  glass  face¬ 
plate,  the  ITO  bias  electrode  and  the  CeOi  blocking  layer  for 
HARP.  The  value  of  Ts  is  a  function  of  the  wavelength  X  and 
was  experimentally  measured  using  a  HARP  substrate  con¬ 
sisting  just  of  the  faceplate,  electrode,  and  blocking  layer. 
During  the  manufacture  of  HARP  targets  this  structure  is 
assembled  prior  to  deposition  of  the  a-Se  layer.  The  substrate 
sample  was  placed  in  the  beam  path  of  the  monochromator 
and  the  transmission  of  light  was  calculated  by  taking  the 
ratio  of  the  light  intensity  measured  by  the  Si  PD  with  the 
sample  in  place  and  that  without. 


Fig.  4.  Optical  quantum  efficiency  of  a- Se  calculated  using  Onsager  theory 
and  r0=  1.7  nm,  compared  with  the  measured  effective  quantum  efficiency 
77*  at  A. =540  nm. 


lanche  multiplication.  To  choose  the  proper  operational  con¬ 
dition  for  HARP  in  different  x-ray  imaging  applications,  we 
need  to  understand  separately  the  ESe  dependence  of  77  and 
gav.  This  will  be  discussed  in  Sec.  II  B. 


B.  Determination  of  quantum  efficiency  and 
avalanche  gain 

The  effective  quantum  efficiency  rf  of  HARP  is  propor¬ 
tional  to  both  77  and  gav.  The  value  of  77  in  a-Se  depends  on 
both  ESe  and  the  wavelength  X  of  the  incident  photon.  It  was 
found  by  Pai  and  Enck  that  the  ESe  and  X  dependence  of 

25 

photogeneration  in  a- Se  follow  the  Onsager  mechanism. 
Onsager  theory  states  that  every  absorbed  photon  creates  a 
pair  of  thermalized  carriers  bound  by  their  Coulombic  attrac¬ 
tion.  The  initial  separation  r0  of  the  charge  pair  increases 
with  incident  photon  energy.  The  probability  of  dissociation 
for  the  bound  pairs  into  free  carriers  increases  with  r0  and  the 
applied  external  field  ESe  with  the  relationship 
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2.  Measurement  results 

The  measured  effective  quantum  efficiency  77  is  shown  in 
Fig.  3  as  a  function  of  ESe  for  an  8  /im  thick  HARP  layer. 
For  blue  light  (X=400nm),  77  reaches  a  plateau  at  ESe 
> 20  V / /ini.  This  corresponds  to  a  saturation  of  optical 
quantum  efficiency  77  at  unity.  With  a  further  increase  in  ESe 
to  >  80  V  /  77, m ,  77*  starts  to  increase  again,  which  shows  the 
onset  of  avalanche.  However  for  green  light,  the  ESe  depen¬ 
dence  of  77  is  quite  different.  As  shown  in  Fig.  3  for  X 
=  540  nm,  which  is  the  peak  of  the  emission  spectrum  of  Csl 
(Tl),  77  starts  at  a  much  smaller  value  compared  to  that  for 
blue  light.  It  then  increases  continuously  with  ESe  without 
reaching  a  plateau  before  avalanche  begins  at  —80  V / /mi. 
This  is  because  the  optical  77  for  green  light  is  lower  and 
does  not  reach  unity  until  ESe  is  well  into  the  range  of  ava- 


with 
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where  tf>0  is  defined  as  the  efficiency  of  production  of  ther¬ 
malized  charge  pairs  per  absorbed  photon  and  is  independent 
of  ESe.  Pai  and  Enck  showed  that  the  room  temperature  mea¬ 
surement  of  r]  can  be  best  fitted  with  the  expression  in  Eq. 
(2)  using  4>{]  =  1  and  r0  in  the  range  between  0.8  and  7.0  nm 
for  X  values  between  620  and  400  nm.  For  the  rest  of  the 
discussion  in  this  paper,  we  will  focus  on  X  =  540  nm,  which 
is  the  peak  of  the  emission  spectrum  for  Csl  (Tl).  Shown  in 
Fig.  4  are  the  values  of  77  calculated  using  Eq.  (2)  and  r0 
=  1.7  nm.  Also  plotted  in  Fig.  4  are  the  77  measured  at  the 
same  wavelength  for  z/Se= 8  /im.  These  results  show  that  the 
Onsager  theory  provides  a  reasonable  prediction  for  the  ESe 
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E^V/pm) 
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1/ESe(Mm/V) 
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Fig.  5.  (a)  Solid  circles  are  the  gav  for  the  8  /im  HARP  layer  calculated  by 
dividing  the  measured  if  at  540  nm  by  the  optical  77  predicted  by  the 
Onsager  theory  using  r0=1.7  nm.  The  solid  line  shows  the  gav  calculated 
using  the  fitted  f}x  and  f32  values,  (b)  The  plot  of  In  [ln(gm)]  as  a  function  of 
1  /ESe,  and  the  best  linear  fit  for  the  data. 


dependence  of  77  in  HARP  for  ESe^10  where  gav=  1.  The 
calculated  77  using  Onsager  theory  was  then  divided  into  the 
measured  77  so  that  the  ESe  dependence  of  avalanche  gain 
gm  can  be  obtained.  The  results  are  shown  in  Fig.  5(a),  where 
gav=46  at  £Sc=  1  10  V/ 77m.  Figure  5(a)  shows  that  gav  in¬ 
creases  very  rapidly  with  ZsSe  hence  choosing  the  appropriate 
ESe  for  HARP  is  crucial.  To  determine  the  x-ray  response  of 
SHARP-AMFPI  it  is  necessary  to  understand  the  continuous 
dependence  of  gdv  on  ESe  because  as  image  charge  accumu¬ 
lates  on  the  pixels,  the  voltage  drop  across  the  HARP  layer 
will  decrease.  This  change  in  ESe  will  produce  a  correspond¬ 
ing  change  in  g.dv  and  affect  the  x-ray  sensitivity,  linearity, 
and  dynamic  range  of  the  FPI. 

The  avalanche  gain  characteristics  in  HARP  have  been 
investigated  as  a  function  of  dSe  and  ESe.26  According  to  the 
theory  of  avalanche  multiplication,  when  ESe  across  the 
thickness  dSe  is  uniform,  and  the  distance  through  which  ava¬ 
lanche  occurs  is  equal  to  dSe  (which  is  true  for  an  8  /im 


HARP  because  photons  are  absorbed  within  0.1  fim  at  the 
entrance  side  of  the  a-Se  layer),  gm  is  given  by"8 

{fi-a)ed 

■?av_  (3  -  aed^iJ3~a)  ’  3 

where  a  and  fi  are  the  impact  ionization  coefficients  of  elec¬ 
trons  and  holes,  respectively.  In  a-Se  the  value  of  fi  is  much 
greater  than  a  since  the  drift  mobility  of  holes 
(0.18  cm2/Vs)  is  much  higher  than  that  of  electrons 
(0.003  cnr/Vs).  This  is  why  in  HARP  light  always  enters 
from  the  positive  bias  electrode  to  maximize  avalanche  gain. 
It  is  important  to  note  that  in  order  for  gdv  [as  defined  in  Eq. 

(3) ]  to  have  a  finite  positive  value,  i.e.,  stable  avalanche,  the 
condition  of  fi>  ae^~°^dse  must  be  satisfied.  Under  the  as¬ 
sumption  that  the  contribution  of  electrons  to  avalanche  gain 
is  negligible,  Eq.  (3)  can  be  simplified  to 

gm  =  ed^  (4) 

with  fi  given  by 

fi=fiie-WE se.  (5) 

The  values  for  fix  and  fi2  can  be  determined  by  curve 
fitting  of  experimental  measurement  of  gdv.  By  taking  the 
logarithms  of  both  sides  of  Eq.  (4)  twice,  we  obtain 

ln(ln  gav)  =  M/Mse)  -  (6) 

ESe 

which  indicates  that  the  quantity  ln(ln  gav)  is  linearly  related 
to  l/ESe.  Shown  in  Fig.  5(b)  is  the  experimentally  deter¬ 
mined  gdv  data  from  Fig.  5(a)  replotted  in  the  form  of 
Inf  In  gav)  as  a  function  of  l/ESe.  The  best  linear  fit  to  the 
experimental  data  is  also  shown  in  Fig.  5(b).  From  the  inter¬ 
cept  and  the  slope  of  the  fitted  line  we  obtained  fix  =  5.5 
X  103  /u,m_1  and  fi2  =  1.029X  103  V//rm.  Using  these  fil  and 
fi2  values  and  the  analytical  expression  for  g.lv  shown  in  Eqs. 

(4)  and  (5),  the  fitted  gdv  values  were  calculated  and  are 
shown  in  Fig.  5(a).  It  shows  that  the  calculated  gav  using  the 
fitted  parameters  agrees  well  with  experimentally  determined 
gm  for  ESe 1 10  V/ /im.  Hence  the  fitted  parameters  fi\  and 
fi2  will  be  used  in  the  following  investigation  of  the  potential 
x-ray  imaging  performance  of  SHARP-AMFPI. 

III.  POTENTIAL  X-RAY  IMAGING  PERFORMANCE 

A.  Cascaded  linear  system  model  for  SHARP- 
AMFPI 

Previously  we  have  developed  a  cascaded  linear  system 

23 

model  for  indirect  FPI  with  avalanche  gain.  This  prelimi¬ 
nary  investigation  did  not  include  K-fluorescence  reabsorp¬ 
tion  in  Csl  or  direct  x-ray  interaction  in  HARP.  We  showed 
that  the  DQE(/)  is  related  to  the  avalanche  gain  gav  through 
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Fig.  6.  Flow  diagram  showing  the  serial  and  parallel 
stages  of  the  cascaded  linear  system  model  for 
SHARP-AMFPI. 


DQE  (/)  = 


T2b{f)T2c(f)T2a{f) 


.  vA 


-TliDtif)  + 


(2gav  ~  1) 

gaVg  Vx  yv\ 


T%f)  + 


APgavgyxyy)  qo 


(7) 


where  q0  is  the  number  of  incident  x-ray  photons  per  unit 
area,  y  is  the  optical  coupling  efficiency  between  Csl  and  the 
HARP  layer,  Th(f),  Tc(f),  and  Tjf)  are  the  MTFs  due  to  blur 
in  Csl,  the  optical  coupling  between  Csl  and  HARP,  and  the 
pixel  aperture  function,  respectively.  As ,  rjx,  and  g  are  the 
Swank  factor,  x-ray  quantum  absorption  efficiency,  and  x-ray 
to  photon  conversion  gain  of  the  Csl  layer,  respectively.  The 
factors  related  to  the  TFT  array  are  the  electronic  noise  Sn 
and  the  pixel  width  ap.  The  denominator  of  Eq.  (7)  is  the 
dose  normalized  NPS.  The  first  term  is  the  x-ray  quantum 
noise,  the  second  term  is  the  secondary  quantum  noise  asso¬ 
ciated  with  the  variance  in  gain  in  the  conversion  and  the 
avalanche  process,  and  the  third  term  is  the  contribution 
from  the  electronic  noise  of  the  TFT  array.  This  expression 
shows  that  the  electronic  noise  is  more  important  at  low 
dose,  and  is  inversely  proportional  to  gav  .  Therefore  intro¬ 
ducing  avalanche  gain  is  an  effective  method  for  reducing 
noise  at  low  dose.  The  secondary  quantum  noise  for  gav>  1 
is  approximately  twice  that  for  gav=l,  assuming  rj  is  the 
same  in  both  cases.  This  excess  secondary  quantum  noise  is 
the  penalty  for  avalanche  gain.  However  for  SHARP-AMFPI 
with  high  conversion  gain  g  and  efficient  optical  coupling 
between  Csl  and  HARP,  i.e.,  grjxyr]>2,  the  secondary  quan¬ 
tum  noise  is  negligible  compared  to  the  x-ray  quantum  noise. 
Hence  the  DQE  of  SHARP-AMFPI  will  not  be  degraded  by 
the  noise  associated  with  avalanche  gain.  This  has  been  dem¬ 
onstrated  experimentally  by  coupling  a  HARP  tube  to  the 
output  of  an  x-ray  image  intensifier  (XRII)  and  measuring 
the  DQE  as  a  function  of  avalanche  gain  gav.  It  was  shown 
that  DQE  remains  essentially  constant  as  a  function  of  gav 
(with  electronic  noise  negligible  in  all  cases). 

In  a  later  investigation  of  the  inherent  imaging  perfor¬ 
mance  of  Csl  layers,  we  included  the  effect  of 
K-fluorescence  reabsorption  of  Csl  in  the  cascaded  linear 
system  model  for  indirect  AMFPI.  Our  results  showed  that 
the  K-fluorescence  reabsorption  in  Csl  has  two  effects:  (1) 
degradation  of  the  Swank  factor  for  x-ray  energies  above  the 
K  edges  and  (2)  additional  image  blur,  which  accounts  for  a 


drop  in  presampling  MTF  of  20-30  %  at  10  cycles /mm  for 
Csl  thickness  between  150  and  600  /l m.  For  the  present  pa¬ 
per,  a  more  complete  cascaded  linear  system  model  was  de¬ 
veloped  that  includes  the  effect  of  direct  x-ray  interaction  in 
HARP.  Shown  in  Fig.  6  are  the  serial  and  parallel  cascaded 
stages  involved  in  the  formation  of  an  x-ray  image  by 
SHARP-AMFPI.  Each  stage  can  be  categorized  as  one  of  the 
following  six  processes:  (1)  gain  or  selection;  (2)  stochastic 
blurring;  (3)  deterministic  blurring;  (4)  aliasing;  (5)  addition; 
and  (6)  parallel  processes.  There  are  two  parallel  processes: 
one  associated  with  the  probability  of  K-fluorescence  genera¬ 
tion  and  reabsorption  in  Csl,  and  the  other  with  direct  inter¬ 
action  of  x  rays  in  the  HARP  layer.  We  have  previously 
described  in  detail  the  transfer  of  spatial  frequency  depen¬ 
dent  signal  '!>(/)  and  NPS,  S(f)  through  all  the  stages  except 
direct  x-ray  interaction  in  the  HARP  layer.  This  process  will 
be  described  in  Sec.  Ill  C  of  this  paper. 

Using  the  cascaded  linear  system  model  shown  in  Figure 
6,  we  investigated  two  important  aspects  of  the  x-ray  imag¬ 
ing  performance  of  SHARP-AMFPI.  One  is  the  pixel  x-ray 
response  of  the  detector,  i.e.,  the  linearity  and  dynamic 
range,  for  both  the  avalanche  and  nonavalanche  operational 
modes;  the  other  is  the  effect  of  direct  x-ray  interaction  in 
HARP  on  the  DQE  of  the  detector.  These  will  be  described 
in  Sections  B  and  C,  respectively. 


B.  Pixel  x-ray  response 

In  this  section,  we  will  apply  the  ESe  dependence  of  gav 
and  rj  derived  previously  in  Sec.  II  B  to  the  cascaded  linear 
system  model  to  determine  the  x-ray  response  of  SHARP- 
AMFPI.  Due  to  the  difference  in  resolution  and  avalanche 
gain  requirement  for  mammography  and  radiography/ 
fluoroscopy  (R/F),  the  calculation  will  be  performed  sepa¬ 
rately  for  these  two  imaging  applications.  From  the  cascaded 
linear  system  model  the  mean  signal  from  each  pixel  of  the 
detector  is  given  by 
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Table  I.  SHARP-AMFPI  detector  design  parameters  and  operating  conditions  chosen  for  different  x-ray  im¬ 
aging  applications. 


Mammography  detector 

R/F  detector 

Tomosynthesis 

Screening 

Fluoroscopy 

Radiography 

Pixel  pitch  ap(/i m) 

70 

200 

X-ray  spectrum 

28  kVp  Mo/Mo,  4  cm  PMMA 

RQA5 

Electronic  noise 

2000  electrons 

1500  electrons 

Csl  (75%  packing) 

150  /xm 

HR 

600  /xm  HL 

W  for  Csl  (eV) 

30 

18 

y 

0.8 

0.8 

Minimum  exposure 

0.1  mR 

1  mR 

0.1  yllR 

30  fiR 

Eseiyu ™) 

110 

105 

110 

95 

V 

0.36 

0.35 

0.36 

0.31 

Sav 

46 

12 

46 

2.3 

apVVgwKb 

w 


(8) 


where  W  is  the  energy  required  to  generate  an  optical  photon 
in  Csl,  and  Eab  is  the  x-ray  energy  absorbed  by  Csl  per  unit 
area.  Eab  includes  contributions  from  all  three  parallel  paths 
associated  with  K-fluorescence 


Eab  =  2  q0(E)  %(£)[(  1  -  Pk)E  +(E- Ek)Pk  +  EkPkfk] , 

E 

(9) 

where  q0(E )  is  the  incident  x-ray  photons  per  unit  area  as  a 
function  of  energy  E  for  a  given  x-ray  spectrum,  Pk  is  the 
probability  of  K  fluorescence,  Ek  is  the  energy  of  the 
K-fluorescence  photon,  and  fk  is  the  fraction  of  K  fluores¬ 
cence  reabsorbed  in  the  Csl  layer.  Each  of  the  K-fluorescence 
related  terms  in  Eq.  (9)  was  determined  separately  for  the  K„ 
and  Kfj  photons  from  both  Cs  and  I  atoms  and  integrated  to 
obtain  the  final  result  of  Eab- 

We  have  previously  investigated  the  inherent  x-ray  imag¬ 
ing  performance  of  Csl  (Tl)  using  individual  samples  opti¬ 
mized  for  either  high  resolution  (HR)  or  high  light  output 
(HL).30  Our  results  showed  that  the  HR  type  has  essentially 
no  depth  dependent  blur  and  as  a  result  provides  better 
DQE(/)  at  high  spatial  frequencies  compared  to  HL  Csl  lay¬ 
ers.  This  was  achieved  at  the  cost  of  a  lower  light  output 
(60%)  for  the  HR  compared  to  the  HL  Csl.  For  the  calcula¬ 
tion  in  this  paper,  we  used  150  pm  HR  Csl  for  mammogra¬ 
phy  and  600  p m  HL  Csl  for  R/F  applications.  Our  previ¬ 
ously  measured  values  for  A  s ,  W,  and  presampling  MTF  of 
these  Csl  layers  were  used  in  the  cascaded  linear  system 
model.  Since  we  have  not  yet  implemented  the  depth  depen¬ 
dent  blur  of  Csl  in  our  cascaded  linear  system  model,  the 
DQE(/)  at  high  frequencies  for  the  600  pm  HL  Csl  layer 
(for  R/F)  will  be  overestimated.  Table  I  summarizes  the  de¬ 
tector  parameters  and  HARP  operating  conditions  chosen  for 
each  application.  For  mammography,  we  consider  two  differ¬ 
ent  image  acquisition  modes:  screening  mammography  and 
tomosynthesis  [where  several  images  of  the  breast  are  ac¬ 
quired  from  different  viewing  angles  for  tomographic  recon¬ 
struction  of  three-dimensional  (3D)  images]. 


1.  Radiography/Fluoroscopy 

Shown  in  Fig.  7(a)  is  the  DQE(/)  calculated  using  the 
cascaded  linear  system  model  for  the  detector  with  different 
gav  settings.  For  comparison  we  also  calculated  the  DQE(/) 


(a) 


(b) 

Fig.  7.  Calculated  DQE  (/)  using  the  detector  parameters  and  operating 
conditions  shown  in  Table  I  for  a  fluoroscopy  detector  at  the  x-ray  exposure 
level  of  (a)  0.1  /xR  and  (b)  30  /xR. 
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(a) 


(b) 

Fig.  8.  Calculated  image  charge  on  each  pixel  of  the  detector  and  the  cor¬ 
responding  avalanche  gain  as  a  function  of  x-ray  exposure  for  the  detector 
operational  setting  of  (a)  fluoroscopy  and  (b)  radiography. 


of  existing  indirect  AMFPI  technology  without  avalanche 
gain  using  77=0.8  [which  is  the  quantum  efficiency  of  a-Si 
PD  (Ref.  1)]  and  an  optimistic  fill  factor  of  unity.  The  calcu¬ 
lation  was  performed  for  the  lowest  exposure  in  fluoroscopy 
(0.1  /rR),  where  the  detector  performance  is  most  suscep¬ 
tible  to  electronic  noise.  It  shows  that  gav=46  is  sufficient  to 
bring  the  DQE(/)  to  its  theoretical  limit  (i.e.,  a  detector  with¬ 
out  electronic  noise).  The  DQE(/)  for  the  next  highest  gav 
=  12  measured  at  ESe=  105  V / /im  also  provides  reasonable 
DQE,  which  is  an  significant  improvement  from  existing  FPI 
technology.  Thus  an  Ese  between  105  and  110  W/rm  will  be 
adequate  for  fluoroscopy.  Figure  7(b)  shows  the  DQE(/)  cal¬ 
culated  for  the  lowest  radiographic  exposure  of  30  /jlR.  It 
shows  that  ESe= 95  V / /xm  can  provide  sufficient  77  (0.31) 
and  gav  (2.3).  Hence  the  ESe  values  chosen  in  Table  I  are 
adequate  for  the  dual-mode  operation  of  a  R/F  detector. 

Shown  in  Fig.  8  are  the  calculated  pixel  x-ray  responses 
for  the  R/F  detector  under  two  different  operational  condi¬ 
tions.  Here  we  assumed  that  the  dominant  pixel  capacitance 
is  the  capacitance  of  the  HARP  layer.  This  condition  results 
in  the  largest  drop  in  the  potential  across  the  HARP  layer  for 
a  given  x-ray  exposure,  and  creates  the  biggest  effect  on 


linearity.  Figure  8(a)  shows  that  for  fluoroscopy,  the  initial 
£’se=H0  V / /im  provides  a  gav  of  46.  Within  the  regular  ex¬ 
posure  range  of  0.1-10  77JR  per  frame  in  fluoroscopy,  the 
detector  has  good  linearity,  i.e.,  an  essentially  constant  gav. 
With  further  increase  in  exposure,  which  corresponds  to  the 
exposure  near  or  beyond  the  edge  of  the  human  body,  the 
x-ray  response  becomes  sublinear  with  a  decrease  in  gm. 
This  decrease  in  gdv  ensures  a  wide  dynamic  range  for  the 
detector  because  there  is  no  “hard”  saturation.  The  pixel  po¬ 
tential  will  continue  to  increase  as  a  function  of  x-ray  expo¬ 
sure  until  it  reaches  the  predetermined  Vmax  for  high  voltage 
protection  of  the  TFTs.  The  method  for  high  voltage  protec- 
tion  was  described  in  detail  elsewhere.  The  pixel  potential 
for  the  exposure  of  1  mR  in  the  fluoroscopy  mode  is  VP 
=  56  V,  which  is  safe  for  properly  designed  a-Si  TFTs.32  The 
mechanisms  for  high  voltage  damage  for  a-Si  TFTs  are  most 
likely  due  to  gate  dielectric  breakdown  or  heating.  In 
AMFPI,  the  readout  of  image  charge  usually  results  in  an 
exponential  decay  of  the  TFT  channel  current,  and  is  not 
likely  to  generate  enough  heating  to  cause  damage.  Hence 
the  main  mechanism  for  high  voltage  damage  in  AMFPI  is 
gate  dielectric  breakdown.  The  breakdown  voltage  depends 
on  the  thickness  and  material  properties  of  the  gate  insulator, 
which  is  usually  silicon  nitride  with  a  thickness  of  several 
thousand  angstroms.34  Since  silicon  nitride  films  can  stand 
an  electric  field  of  up  to  103  V  / /im.  a  pixel  potential  of  up  to 
100  V  may  be  safe  for  the  TFTs.  Figure  8(b)  shows  that  with 
a  much  lower  gav=2.3  programmed  for  radiography,  the  de¬ 
tector  response  for  the  regular  exposure  range  of 
30  /t.R-3  mR  is  essentially  linear  with  very  little  change  in 
gav.  Comparing  Figures  8(a)  and  8(b),  the  image  charge  on 
the  pixel  is  in  the  same  range  despite  their  large  difference  in 
exposure  range.  This  demonstrates  the  wide  dynamic  range 
(over  five  orders  of  magnitude)  of  the  detector  that  is  made 
possible  by  varying  gav  for  different  x-ray  imaging  applica¬ 
tions.  This  is  an  advantage  over  detectors  with  constant  high 
gain,  which  have  the  problem  of  pixel  saturation  for  radio- 
graphic  applications. 

2.  Mammography 

Figures  9(a)  and  9(b)  show  the  DQE  calculated  for  mam¬ 
mography  using  the  detector  parameters  and  operating  con¬ 
ditions  listed  in  Table  I  for  tomosynthesis  and  screening  ap¬ 
plications.  Also  shown  in  Fig.  9(a)  is  the  calculated  DQE  of 
indirect  FPI  without  avalanche  gain  (77=0.8  and  unity  fill 
factor),  where  the  electronic  noise  causes  severe  degradation 
of  DQE  at  the  low  dose  for  each  view  of  tomosynthesis. 
Figure  9  shows  that  the  choices  of  ESe  for  both  tomosynthe¬ 
sis  and  screening  mammography  modes  are  adequate  to 
achieve  the  theoretical  limit  of  DQE  without  electronic 
noise. 

The  pixel  x-ray  response  for  mammography  was  calcu¬ 
lated  using  Eq.  (8)  and  the  results  are  shown  in  Figs.  10(a) 
and  10(b)  for  tomosynthesis  and  screening  mammography, 
respectively.  Two  sets  of  curves  are  shown  in  each  graph  of 
Fig.  10,  one  without  additional  pixel  storage  capacitance,  and 
the  other  with  a  storage  capacitor  of  Cst=0.5  pF.  The  mini- 
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Spatial  frequency  (cycles/mm) 
(a) 


Spatial  frequency  (cycies/mm) 

(b) 

Fig.  9.  Calculated  DQE  (/)  using  the  detector  parameters  and  operating 
conditions  shown  in  Table  I  for  mammography  (a)  tomosynthesis  with  mini¬ 
mum  exposure  of  0.1  mR  and  (b)  screening  with  minimum  exposure  of  1 
mR. 


mum  exposure  in  each  graph  of  Fig.  10  corresponds  to  the 
detector  exposure  behind  dense  breast  tissue,  and  the  maxi¬ 
mum  corresponds  to  raw  exposure.  Figure  10  shows  that 
when  there  is  no  additional  storage  capacitance,  the  signal 
increases  sublinearly  as  a  function  of  the  x-ray  exposure. 
This  is  because  the  pixel  capacitance  due  to  the  HARP  layer 
for  mammography  is  only  —0.03  pF  with  cip=lQ  /xm,  which 
is  much  smaller  than  that  for  the  R/F  detector  (0.28  pF). 
Hence  the  signal  accumulated  on  each  pixel  creates  a  much 
more  significant  drop  in  ESe,  which  causes  a  rapid  decrease 
in  gav.  The  maximum  image  charge  shown  in  Figs.  10(a)  and 
10(b)  corresponds  to  a  pixel  potential  of  —100  and  210  V, 
respectively.  These  potentials  exceed  the  maximum  safe  volt¬ 
age  for  the  TFTs.  One  method  for  alleviating  this  problem 
while  maintaining  the  same  dynamic  range  is  to  add  a  pixel 
storage  capacitance  Cst.  With  Cst=0.5  pF,  the  pixel  potential 
rises  much  more  slowly  and  the  x-ray  response  of  the  detec¬ 
tor  becomes  more  linear,  as  shown  in  Fig.  10.  The  pixel 


(a) 


Fig.  10.  Calculated  image  charge  on  each  pixel  of  the  mammography  de¬ 
tector  and  the  corresponding  avalanche  gain  as  a  function  of  x-ray  exposure 
with  or  without  the  addition  of  a  pixel  storage  capacitor.  The  results  are  for 
detector  operating  conditions  chosen  for  (a)  tomosynthesis  image  acquisi¬ 
tion  and  (b)  screening  mammography. 


potential  corresponding  to  the  maximum  exposures  in  Figs. 
10(a)  and  10(b)  is  only  25  and  50  V,  respectively,  both  of 
which  are  safe  for  the  TFTs. 

C.  Direct  x-ray  interaction  in  HARP 

In  indirect  AMFPI  using  phosphor  and  a-Si  photodiodes, 
direct  interaction  of  x  rays  in  the  photodiodes  produces  ad- 
ditional  noise.  This  is  because  the  x  ray  to  charge  conver¬ 
sion  gain  in  the  photodiodes  is  much  higher  (e.g.,  45  times  in 
the  example  of  Ref.  35)  than  that  in  the  phosphors,  and  there¬ 
fore  each  direct  interaction  x  ray  produces  much  higher  sig¬ 
nal  (and  noise)  than  that  absorbed  in  the  phosphor.  For  our 
proposed  detector,  the  HARP  sensor  (e.g.,  t/Se=8  /j.iti)  is 
much  thicker  than  the  a-Si  photodiodes  (—  1  /i m)  in  existing 
indirect  FPI,  hence  there  will  be  more  direct  x-ray  interaction 
events.  Fortunately  the  x  ray  to  charge  conversion  gain  in 
HARP36  ( VFSe  =  1 6  eV  at  £Se=100  V)  is  very  similar  to  that 
in  Csl,  thus  direct  x-ray  interaction  in  HARP  is  not  expected 
to  add  additional  noise  at  low  spatial  frequencies.  At  high 
spatial  frequencies,  however,  the  effect  of  direct  x-ray  inter¬ 
action  could  be  more  significant.  This  is  because  both  the 
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Indirect  Direct 


Fig.  1 1 .  Diagram  showing  the  depth  dependent  gain  for  direct  x-ray  inter¬ 
action  in  HARP.  For  the  x  rays  absorbed  in  Csl,  the  light  photons  created 
enter  the  top  interface  of  HARP,  where  holes  are  generated  and  undergo 
avalanche  multiplication  for  the  entire  thickness  of  HARP  dSe.  For  the  x  rays 
interacting  directly  with  HARP  at  depth  x  from  the  top  interface,  the  holes 
generated  undergo  avalanche  for  the  thickness  of  (dSe-x),  which  causes 
depth  dependence  avalanche  gain  variation. 


signal  and  the  NPS  due  to  x  rays  absorbed  in  Csl  decrease  as 
a  function  of  spatial  frequency  due  to  blur  in  Csl,  whereas 
the  signal  and  NPS  due  to  direct  x-ray  interaction  in  HARP 
are  essentially  white  (without  including  the  effect  of  K  fluo¬ 
rescence).  Furthermore,  when  HARP  is  operated  in  the  ava¬ 
lanche  mode,  the  x  rays  absorbed  by  HARP  will  experience 
different  avalanche  gain  depending  on  their  depth  of  interac¬ 
tion.  As  shown  in  Fig.  1 1 ,  x  rays  absorbed  in  HARP  at  a 
depth  x  from  the  top  (light  entrance)  interface  have  an  effec¬ 
tive  thickness  for  avalanche  multiplication  of  (c/Se-x).  This 
could  lead  to  a  depth  dependent  variation  in  gav,  and  add 
additional  noise/  In  the  following  we  will  derive  the  propa¬ 
gation  of  signal  and  noise  due  to  direct  x-ray  interaction  in 
HARP  using  the  How  diagram  shown  in  Fig.  6,  and  investi¬ 
gate  its  effect  on  the  DQE  of  SHARP-AMFPI. 

1.  Signal  and  noise  propagation  for  direct  x-ray 
interaction 

As  shown  in  Fig.  6,  the  input  to  the  parallel  path  for  direct 
x-ray  interaction  in  HARP  is  the  x  rays  transmitted  through 
the  Csl  layer  q h  which  is  a  function  of  x-ray  energy.  For 
clarity  we  assume  a  monoenergetic  x-ray  beam  in  the  follow¬ 
ing  derivation  of  theory,  although  we  used  the  full  x-ray 
spectrum  in  our  calculation  of  direct  x-ray  interaction.  We 
ignored  the  effect  of  K-fluorescence  reabsorption  in  HARP. 
Under  this  assumption,  both  the  signal  and  NPS  for  direct 
x-ray  interaction  are  white  prior  to  the  pixel  aperture  func¬ 
tion,  hence  we  could  omit  the  frequency  dependence  in  the 
following  derivation.  The  x-ray  signal  produced  at  depth  x 
into  the  HARP  layer  by  a  small  thickness  dx  of  a-.Se  is  given 
by 

$(*)  =  <7 1  e-/“gav(x)  gsePdx ,  (10) 

where  p  is  the  linear  attenuation  coefficient  of  a-Se,  and  gSe 
is  the  average  x  ray  to  charge  conversion  gain  which  is  equal 


to  £7lTse-  The  depth  dependence  in  gav  can  be  derived  from 
Eq.  (4)  as  gav(x)  =exp[/3((fSe-x)].  The  NPS  due  to  direct  in¬ 
teraction  as  a  function  of  depth  S(x)  is  given  by37 

S (x)  =  q  i  e-/ir[  (g2Se  +  gavM  +  gs^Jx)  ]  p  dx .  (11) 

Integrating  both  Eqs.  (10)  and  (11)  with  respect  to  x  from  0 
to  dSe,  we  obtain  the  total  signal  and  NPS  due  to  direct  x-ray 
interaction 

^_qxpgsJ,e^-e-^) 

p  +  /3 


and 

q\  Msie  +  °L)(e2Mie " 

p  +  2f3 


(13) 


2.  Calculation  results 

Using  Eqs.  (12)  and  (13),  we  calculated  the  contribution 
to  the  total  signal  and  NPS  by  the  direct  interaction  of  x  rays 
in  HARP  using  the  parameters  shown  in  Table  I  for  a  R/F 
detector.  With  a  RQA5  x-ray  spectrum  and  a  600  pm  HL 
layer  of  Csl,  the  fraction  of  x  rays  transmitted  through  the 
Csl  layer  is  13.2%,  out  of  which  —1%  is  absorbed  by  the 
HARP  layer.  Shown  in  Fig.  12(a)  is  a  comparison  between 
the  presampling  signal  spectra  (before  the  pixel  aperture 
function)  due  to  x  rays  absorbed  in  Csl  and  direct  x-ray 
interaction  in  HARP.  At  low  frequencies,  the  signal  due  to 
direct  x-ray  interaction  is  less  than  0.2%  of  that  due  to  x-ray 
absorbed  in  Csl.  This  is  consistent  with  the  difference  in 
x-ray  quantum  efficiency  between  Csl  and  HARP.  At  high 
frequencies  (e.g.,  5  cycles/mm),  the  signal  from  Csl  de¬ 
creases  due  to  the  blur  in  Csl,  however  the  signal  from  direct 
x-ray  interaction  in  HARP  stays  constant.  But  the  latter  is 
still  negligible.  Shown  in  Fig.  12(b)  is  the  NPS  (before  pixel 
aperture  function)  comparison.  At  5  cycles/mm,  the  NPS 
due  to  direct  x-ray  interaction  becomes  significant  (—17%  of 
the  NPS  due  to  Csl).  However  below  the  Nyquist  frequency 
of  2.5  cycles/mm,  it  is  still  negligible  (<3%).  This  is  why 
there  is  no  significant  change  in  DQE  after  including  the 
direct  x-ray  interaction  events  in  HARP,  as  shown  in  Fig.  13. 
Therefore  we  can  conclude  that  direct  x-ray  interaction  in 
HARP  in  an  indirect  FPI  with  avalanche  gain  has  no  degra¬ 
dation  effect  on  DQE  of  the  detector. 

IV.  EFFECT  OF  HARP  THICKNESS  UNIFORMITY 

An  important  issue  for  making  a  large  area  HARP  is  its 
thickness  uniformity.  Currently  HARP  films  can  be  made 
very  uniform  for  an  area  of  1-2  in.  in  diameter  with  essen¬ 
tially  no  visible  gain  variation  in  the  images.  The  images 
currently  produced  by  HARP  for  broadcast  applications  do 
not  require  any  gain  uniformity  correction.  The  area  of 
HARP  needs  to  be  increased  significantly  for  x-ray  imaging 
applications.  Uniform,  large  area  a-Se  films  have  already 
been  developed  for  direct  FPI,  however  the  uniformity  re¬ 
quirement  for  HARP  in  indirect  FPI  with  avalanche  gain  is 
expected  to  be  much  higher.  In  this  section  we  will  derive  the 
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(a) 


Spatial  frequency  (cycles/mm) 

(b) 

Fig.  12.  Comparison  of  the  presampling  signal  and  NPS  (before  pixel  ap¬ 
erture  function)  due  to  x  ray  absorbed  in  Csl  and  direct  x-ray  interaction  in 
HARP:  (a)  signal  spectra,  where  the  signal  due  to  direct  interaction  is  neg¬ 
ligible  and  (b)  NPS,  where  the  NPS  due  to  direct  interaction  under  the 
Nyquist  frequency  of  2.5  cycles/mm  is  negligible. 


Fig.  13.  Comparison  between  the  DQE  of  SHARP-AMFPI  with  and  with¬ 
out  consideration  of  direct  x-ray  interaction  in  HARP. 


relationship  between  thickness  nonuniformity  and  gain  non¬ 
uniformity  for  HARP  sensors  so  that  the  results  can  be  used 
as  a  guideline  for  engineering  of  large  area  HARP  layers  for 
FPL 


A.  Relationship  between  HARP  avalanche  gain  and 
thickness  uniformity 

During  the  operation  of  HARP,  a  constant  potential  Vb  is 
applied  to  the  top  ITO  bias  electrode,  as  shown  in  Fig.  1(b). 
After  each  image  readout,  the  TFT  array  resets  the  potential 
on  pixel  electrodes  to  Vp=0  V.  A  nonuniformity  in  the  thick¬ 
ness  of  HARP  will  result  in  a  variation  in  ESe  and  hence  a 
spatial  nonuniformity  in  gain.  An  increase  in  dSe  under  a 
constant  Vh  causes  a  decrease  in  ESe.  Although  on  one  hand 
an  increase  in  dSe  tends  to  increase  gav  [Eq.  (4)],  it  also 
results  in  a  simultaneous  decrease  in  ESe  under  a  constant  Vb. 
The  overall  effect  is  a  decrease  in  gav  because  of  its  addi¬ 
tional  exponential  dependence  on  ESe  shown  in  Eq.  (5).  By 
substituting  ESe=  Vh/d$e  into  Eq.  (6)  and  differentiating  both 
sides,  we  obtain 


1  Agav  AdSe  /32  A  , 

1ft  <?av  <?av  ^Se  ^ b 


Equation  (14)  can  be  rewritten  as 


gav 


:ln(gj  1 


Pi  \ 

Ese!  dSe 


(14) 


(15) 


Equation  (15)  shows  that  with  (32>ESe ,  which  is  true  for 
HARP,  a  positive  change  in  <fSe  results  in  a  negative  change 
in  gav.  Equation  (15)  also  shows  that  the  percentage  variation 
in  avalanche  gain  is  proportional  to  ln(gav)  and  the  percent¬ 
age  thickness  variation,  hence  HARP  layers  operated  with 
higher  gav  are  more  susceptible  to  gain  nonuniformity  caused 
by  thickness  nonuniformity. 


B.  Calculation  of  gain  nonuniformity 

Using  Eqs.  (4)  and  (5),  and  the  values  for  /3t  and  f32 
determined  in  Sec.  II,  the  variations  in  gav  were  calculated 
for  both  positive  and  negative  change  in  dSe  around  8  /xm, 
and  results  are  shown  in  Fig.  14.  Three  Vb  values  of  840, 
864,  and  880  V  were  used,  which  corresponds  to  gav  of  12, 
25,  and  46,  respectively.  Figure  14  shows  that  the  slope  of 
the  curve,  i.e.,  percentage  change  in  gav,  is  much  higher  for 
negative  change  in  c/Se,  which  causes  an  increase  in  gav,  than 
it  is  for  positive.  Hence  when  choosing  the  bias  potential  Vb 
of  HARP  for  a  desired  ESe ,  the  smallest  dSe  for  the  entire 
detector  area  should  be  used,  so  that  the  variation  in  thick¬ 
ness  is  positive  from  the  reference  thickness.  If  we  assume 
the  smallest  afSe  =  8  /ini,  the  percentage  changes  in  gav  due  to 
an  increase  in  dSe  of  1%  are  19,  23,  and  26  %,  respectively, 
for  the  three  Vb  settings  of  840,  864,  and  880  V.  The  change 
in  gdv  increases  to  33,  40,  and  44  %,  respectively,  when  the 
variation  in  dSe  increases  to  2%.  This  range  of  gav  nonunifor¬ 
mity  is  certainly  within  the  capability  of  gain  correction  and 
the  dynamic  range  of  the  electronic  circuits  of  FPL  With  a 
further  increase  in  dSe  to  4%,  the  value  for  gm  will  decrease 
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<4>m) 


Fig.  14.  Avalanche  gain  gau  calculated  as  a  function  of  the  thickness  of 
HARP,  which  varies  around  8  yum,  under  a  constant  bias  potential  Vb  of  840, 
864,  and  880  V. 


from  the  original  setting  of  46,  25,  and  12  to  15,  10,  and  6, 
respectively,  which  is  a  gain  nonuniformity  of  two  to  three 
folds.  Although  the  variation  in  gav  is  significant,  it  is  not 
expected  to  cause  significant  degradation  in  the  DQE  perfor¬ 
mance,  as  can  be  seen  in  the  example  of  gav= 46  and  12 
shown  in  Fig.  7(a).  However  this  will  pose  a  greater  chal¬ 
lenge  for  the  uniformity  correction  algorithm.  Hence  we  can 
conclude  that  a  thickness  nonuniformity  of  2%  is  desirable 
for  a  HARP  layer,  and  4%  nonuniformity  could  still  possibly 
be  tolerated. 

V.  CONCLUSIONS 

The  feasibility  of  a  concept  of  indirect  conversion  flat- 
panel  detector  with  avalanche  gain  is  being  investigated.  The 
avalanche  gain  and  quantum  efficiency  of  an  8  /jl m  thick 
HARP  layer  were  investigated  experimentally  as  function  of 
ESe.  These  characteristics  of  HARP  were  applied  to  a  cas¬ 
caded  linear  system  to  determine  the  optimal  operating  con¬ 
ditions  required  for  x-ray  imaging.  Our  results  showed  that 
by  varying  the  avalanche  gain  between  2  and  46,  x-ray  quan¬ 
tum  noise  limited  performance  and  a  wide  dynamic  range 
can  be  obtained  for  both  R/F  and  mammographic  imaging 
applications.  Direct  x-ray  interaction  in  HARP  was  shown 
not  to  affect  the  imaging  performance.  Due  to  the  rapid  in¬ 
crease  in  gav  as  a  function  of  ESe,  a  small  variation  in  dSe  can 
cause  significant  nonuniformity  in  avalanche  gain.  It  is  nec¬ 
essary  to  keep  the  thickness  uniformity  of  a  HARP  layer 
within  4%,  which  could  maintain  a  high  DQE  across  the 
entire  detector  and  make  the  variation  in  gav  correctable  us¬ 
ing  gain  calibration. 
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Indirect  flat  panel  imagers  have  been  developed  for  digital  radiography,  fluoroscopy  and  mammog¬ 
raphy,  and  are  now  in  clinical  use.  Screens  made  from  columnar  structured  cesium  iodide  (Csl) 
scintillators  doped  with  thallium  have  been  used  extensively  in  these  detectors.  The  purpose  of  this 
article  is  to  investigate  the  effect  of  screen  optics,  e.g.,  light  escape  efficiency  versus  depth,  on  gain 
fluctuation  noise,  expressed  as  the  Swank  factor.  Our  goal  is  to  obtain  results  useful  in  optimizing 
screens  for  digital  radiography  systems.  Experimental  measurements  from  structured  Csl  samples 
were  used  to  derive  their  screen  optics  properties,  and  the  same  methods  can  also  be  applied  to 
powder  screens.  Csl  screens,  all  of  the  same  thickness  but  with  different  optical  designs  and 
manufacturing  techniques,  were  obtained  from  Hamamatsu  Photonics  Corporation.  The  pulse 
height  spectra  (PHS)  of  the  screens  were  measured  at  different  x-ray  energies.  A  theoretical  model 
was  developed  for  the  light  escape  efficiency  and  a  method  for  deriving  light  escape  efficiency 
versus  depth  from  experimental  PHS  measurements  was  implemented  and  applied  to  the  Csl 
screens.  The  results  showed  that  the  light  escape  efficiency  varies  essentially  linearly  as  a  function 
of  depth  in  the  Csl  samples,  and  that  the  magnitude  of  variation  is  relatively  small,  leading  to  a  high 
Swank  factor.  ©  2006  American  Association  of  Physicists  in  Medicine.  [DOI:  10.1118/1.21 88082] 

Key  words:  indirect  detection,  flat-panel  detectors,  digital  radiography,  cesium  iodide,  Swank 
factor 


I.  INTRODUCTION 

Indirect  flat  panel  imagers  (FPI)  have  been  developed  for 
digital  radiography  (DR),  fluoroscopy,  and  mammography, 
and  are  now  in  clinical  use.  These  imagers  use  a  scintillating 
x-ray  screen  as  the  primary  x-ray  detector  and  read  the  opti¬ 
cal  signal  electronically  with  a  large  area  active  matrix  array 
with  an  amorphous  silicon  (a-Si)  photodiode  and  a  switching 
thin  film  transistor  at  each  pixel.  Although  powdered  phos¬ 
phor  screens,  like  those  used  for  many  years  in  conventional 
screen/film  radiography  have  been  employed,1  structured 
(needle-like)  screens  using  cesium  iodide  (Csl)  are  more 
widely  used  in  FPI.  The  Csl  used  in  FPI  are  doped  with 

thallium  (Tl),  which  is  different  from  the  sodium  (Na)  doped 

2 

Csl  used  in  image  intensifies. 

Because  in  DR  the  detection,  processing,  and  display 
functions  are  independent,  they  can  be  optimized  separately. 
The  detective  quantum  efficiency  (DQE)  has  been  used  as 
the  fundamental  image  quality  metric  for  comparison  be¬ 
tween  and  optimization  of  DR  detectors.  Both  theoretical 
modeling  and  experimental  characterization  of  DQE  have 
been  performed  for  both  direct3  and  indirect4  FPI.  Different 
types  and  thicknesses  of  screens  have  been  used  in  indirect 
FPI  for  different  x-ray  imaging  applications.  It  is  important 
to  understand  and  quantify  the  tradeoffs  between  x-ray  ab¬ 
sorption  efficiency  and  resolution  of  different  screen  designs 


and  their  effects  on  DQE.  The  effect  of  the  screen  component 
of  a  DR  detector  on  the  DQE  has  received  attention  recently 
for  both  powder5  and  structured6  screens. 

The  study  of  the  optical  properties  of  x-ray  screens  pre¬ 
dates  their  use  in  modern  digital  systems.  A  classic  work  by 
Swank  quantified  the  effect  of  variability  in  the  size  of  the 
scintillation  pulses  on  DQE.  This  variability  is  called  gain 
fluctuation  noise,  or  “Swank  noise.”  Swank  also  calculated 
the  blurring  effect  of  light  scattering  and  the  light  escape 
efficiency  in  a  turbid  (light  scattering)  x-ray  screen,8  making 
use  of  previous  work  by  Marshak  et  al.,9  who  was  studying 
neutron  diffusion,  and  by  Hamaker,10  who  had  calculated  the 
light  emitted  from  a  layer  of  x-ray  fluorescent  material.  Ha¬ 
maker  had  used  the  well  known  multiple  light  scattering 
theory  of  Kubelka  and  Munk,11  which  in  turn  was  based  on 
the  equations  developed  by  Schuster  for  light  propagation 
in  stellar  atmospheres. 

In  this  article,  we  extend  some  of  the  previous  work  to 
investigate  the  effect  of  screen  optics,  e.g.,  light  escape  effi¬ 
ciency  versus  depth,  on  Swank  factor,  which  is  proportional 
to  DQE  of  the  detector  at  zero  spatial  frequency.  Our  goal  is 
to  obtain  results  useful  in  optimizing  screens  for  DR  sys¬ 
tems.  As  an  example,  experimental  measurements  from 
structured  Csl  phosphors  were  used  to  derive  their  screen 
optics  properties,  but  the  same  methods  can  be  applied  to 
powder  screens. 
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II.  THEORY  AND  BACKGROUND  FOR  SCREEN 
OPTICS 

A.  Escape  efficiency  (e) 

In  a  phosphor  screen,  not  all  light  photons  generated  by  x 
rays  escape  the  screen  to  contribute  to  signal.  As  shown  in 
Fig.  1,  we  assume  that  the  x-ray  photons  enter  through  the 
backing  of  the  screen,  and  an  optical  sensor  is  placed  next  to 
the  front  surface.  We  assume  that  a  number  m  of  photons  are 
created  in  the  phosphor  layer  of  the  screen,  either  by  prompt 
luminescence  after  x-ray  absorption  in  a  DR  screen,  or  by 
stimulated  luminescence  in  a  computed  radiography  (CR) 
screen.  In  a  nontransparent  screen,  these  photons  will  scatter 


from  discontinuities  in  the  optical  index  of  refraction  in  the 
material,  and  may  also  be  absorbed  in  the  bulk  or  the  back¬ 
ing  of  the  screen.  On  average,  n  photons  will  make  their  way 
to  the  front  (exit)  surface,  where  the  optical  sensor  is  placed, 
and  escape.  We  call  the  fraction  n/m  the  escape  efficiency  or 
e.  It  depends  on  the  thickness,  reflectance  of  both  surfaces 
and  optical  properties  of  the  screen,  and  on  the  depth  of 
creation  of  the  original  photons.  Lubinsky  et  al.  developed 
an  explicit  formula  for  the  escape  efficiency,  e,  of  an  isotro¬ 
pic  powder  screen  in  the  context  of  a  systems  model  for  CR, 
using  the  Kubelka-Munk  theory  and  summing  multiple  scat¬ 
tering  terms 


e(z,L,S,K,Rt,Rb) 


1  -  Rt\  I  b(  1  +  Rb)cosh[bS(L  -  z)]  +  (1  -  Rb)(  1  +  a)sinh[i'5'(L  -  z)] 
2  / 1  b(l  -  RtRb)cosh(bSL)  +  [a(1  +  RtRb)  -  (Rt  +  Rb)]sinh(bSL) 


(1) 


where  z  is  the  depth  of  light  generation  with  z  =  0  represent¬ 
ing  the  light  exit  side  of  the  screen,  L  is  the  thickness  of  the 
screen,  S  and  K  are  the  scattering  and  absorption  coefficients 
in  the  Kubelka-Munk  theory,  Rt  and  Rb  are  the  reflectances 
at  the  top  and  bottom  surfaces,  and  a  and  b  are  defined  as 
fl=  1  +K/S,b  =  \lci2-  1,  respectively.  It  should  be  noted  that 
the  optical  coefficients  S  and  K  do  not  bear  a  simple  relation¬ 
ship  to  the  intrinsic  optical  constants  of  the  screen  medium 
or  the  physical  parameters  such  as  the  particle  size,  as  is 
pointed  out  in  a  recent  paper  on  the  interpretation  of  these 
parameters.14  It  is  not  the  intention  of  our  present  work  to 
derive  the  values  of  S  and  K,  but  rather  to  deduce  the  escape 
efficiency  from  experimental  characterization  of  the  screen. 


A 


L 


V 


Back 


X-ray  photon 


Fig.  1.  Illustrating  the  creation  and  escape  of  photons  from  a  scattering 
medium. 


The  escape  efficiency  model  given  in  Eq.  (1)  can  be  used 
to  study  the  effects  of  screen  optics  design.  Figure  2  shows 
examples  of  e,  calculated  as  a  function  of  depth  of  light 
generation  using  Eq.  (1),  for  different  screen  optics  designs. 
We  refer  to  the  graph  of  escape  efficiency  e  versus  depth  z  as 
the  e-z  curve.  In  Fig.  2(a),  the  light  absorption  coefficient  in 
the  bulk  of  the  screen  is  varied,  and  the  e-z  curve  is  plotted 
for  a  highly  scattering  material,  with  perfectly  nonreflecting 
boundary  conditions  at  both  surfaces  of  the  screen  (back  and 
front).  The  e  is  seen  to  vary  strongly  with  depth,  and  as  the 
absorption  is  increased,  the  decrease  in  e  as  a  function  of  z 
becomes  more  rapid.  Thus,  although  adding  bulk  absorption 
to  a  screen  will  limit  the  extent  of  spatial  light  spreading, 
which  improves  the  MTF,  it  will  also  reduce  the  depth  from 
which  optical  signal  can  propagate  to  the  optical  detector, 
thus  reducing  the  noise-effective  thickness  of  the  screen.  In 
Fig.  2(b),  the  reflectance  of  the  screen  backing  is  varied.  It  is 
noted  that  even  a  backing  reflectance  as  high  as  90%,  is  not 
very  effective  in  boosting  e  for  x  rays  absorbed  near  the 
backing.  This  is  because  in  a  highly  scattering  medium,  a 
photon  created  near  the  back  has  typically  a  lot  more  than 
one  chance  to  be  absorbed  as  it  diffuses  near  the  back  reflec¬ 
tor.  In  order  to  appreciably  improve  e,  the  backreflectance 
must  be  near  unity.  Finally,  in  Fig.  2(c),  the  front  (exit)  sur¬ 
face  reflectance  is  varied.  A  nonzero  exit  surface  reflectance 
might  occur  due  to  internal  reflection  at  the  boundary  when 
the  optical  coupling  between  the  exit  surface  and  the  detector 
is  imperfect.  Here  e  isroughly  linear  within  the  bulk  of  the 
screen,  but  lower  than  unity  at  all  depths  due  to  internal 
reflection  and  absorption  losses. 


B.  Swank  factor  and  pulse  height  spectrum 

The  zero  spatial  frequency,  or  large  area  limit  of  the  DQE, 
DQE(0),  for  a  radiographic  screen  is  given  by 
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Fig.  2.  (a)  Escape  efficiency  vs  depth.  Scattering  coefficient=  100/mm.  Top 
and  bottom  reflectance =0.  Absorption  coefficient  varied:  0.001,  0.01,  0.1, 
1.0/mm.  (b)  Escape  efficiency  vs  depth.  Scattering  coefficient  =  100 /mm. 
Absorption  coefficient =0.001 /mm.  Top  reflectance =0.  Bottom  reflectance 
varied:  0,  0.5,  0.9,  0.98,  0.99.  (c)  Escape  efficiency  vs  depth.  Scattering 
coefficient=  100/mm.  Absorption  coefficient= 0.001 /mm.  Bottom 
reflectance =0.95.  Top  reflectance  varied:  0,  0.75,  0.5. 


DQE(0)  =  V As,  (2) 

where  rj  is  the  x-ray  absorption  efficiency  and  As  is  the  scin¬ 
tillation  efficiency  of  the  screen  (which  is  often  referred  to  as 
the  Swank  factor,  and  sometimes  also  denoted  /).  In  an  inte¬ 
grating  detector,  As  arises  because  the  size  of  the  optical 
pulses  produced  by  absorbed  x  rays  is  not  a  delta  function 
but  has  a  distribution.  (In  a  detector  with  photon  counting 
capability,  the  noise  factor  As  due  to  varying  pulse  sizes  is 
absent.)  Swank7  decomposed  the  optical  pulse  distribution 
into  components  due  to  variations  in  optical  propagation 
(OPD),  variations  in  absorbed  energy  (AED),  and  the  energy 
span  of  the  incident  x-ray  spectrum.  As  can  be  expressed  in 
terms  of  moments  of  the  scintillation  spectrum,  i.e.,  the  prob¬ 
ability  distribution  function  of  the  size  of  the  optical  pulses 
resulting  from  each  x-ray  interaction,  as 


7 

Swank  showed  that  for  monoenergetic  x  rays,  and  under 
certain  reasonable  conditions  for  the  optical  pulse  distribu¬ 
tion  (OPD), 

A-s  =  Adpd^aed  •  (4) 

We  refer  to  the  measured  scintillation  spectrum,  using  mo¬ 
noenergetic  incident  x  rays,  as  the  pulse  height  spectrum 
(PHS).  The  PHS  has  been  measured  for  a  variety  of  practical 
x-ray  screens  using  photon  counting,15'16  pulse  shaping  and 
multichannel  pulse  height  analysis,  and  by  synchronous  de- 
tection  of  a  chopped  x-ray  source.  It  has  previously  been 
shown  experimentally  that  As  is  essentially  constant  with  en¬ 
ergy  below  the  K  edge  of  the  screen,  and  drops  to  a  mini¬ 
mum  just  above  the  K  edge,  due  to  the  escape  of  character- 

17 

istic  K  fluorescence. 

Below  the  K  edge,  and  for  monoenergetic  x  rays,  most  of 
the  variation  observed  in  the  measured  PHS  is  due  to  varia¬ 
tions  in  e  from  various  depths  of  creation  in  the  x-ray  screen. 
By  measuring  the  PHS,  we  obtain  experimental  information 
that  can  potentially  be  used  to  estimate  the  variation  of  e 
with  depth  in  the  scintillating  screen  material. 

III.  MATERIAL  AND  METHODS 
A.  Material 

We  have  investigated  the  inherent  imaging  performance 
of  Csl  with  different  thicknesses  and  types  for  optimization 
of  indirect  FPI  with  and  without  avalanche  gain.19  Here  we 
included  the  two  150  /am  thick  screens  from  the  previous 
study  and  added  two  new  screens  of  the  same  thickness  but 
with  entirely  different  manufacturing  technique.  All  samples 
were  manufactured  by,  and  are  available  commercially  from 
Hamamatsu  Photonics  Corporation.  Aside  from  the  differ¬ 
ence  in  substrate  materials  and  treatment  of  the  surfaces,  the 
intrinsic  Csl  material  properties  (e.g.,  doping)  is  the  same  for 
all  samples.  As  shown  in  Fig.  3,  the  previous  two  Csl  layers 
were  deposited  on  fiber  optic  faceplate  (FOP),  with  metallic 
and  protective  coating  deposited  on  top  of  Csl.  These  layers 
were  called  fiber  optic  scintillator  (FOS).  During  x-ray  expo- 
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Fig.  3.  Diagram  illustrating  four  screen  samples. 

sure,  the  protective  coating  side  of  the  Csl  layer  faced  the 
x-ray  source,  and  the  FOP  side  was  coupled  to  an  optical 
sensor,  e.g.,  a  photomultiplier  tube  (PMT).  Two  types  of 
150  yum  FOS  were  used  in  our  experiment,  one  optimized 
for  high  light  (HL)  output,  and  the  other  optimized  for  high 
resolution  (HR).  Both  HR  and  HL  layers  had  metallic  back¬ 
ing,  but  the  HR  type  had  different  bulk  material  properties  to 
minimize  light  spreading  between  columns,  which  improved 
resolution  at  the  cost  of  lower  light  output  (60%  of  that  for 
HL).  The  two  new  samples  added  for  the  present  study  had 
the  Csl  deposited  on  amorphous  carbon  (a-C)  substrate, 
hence,  called  amorphous  carbon  scintillator  (ACS).  They 
also  came  with  HR  and  HL  types.  Compared  to  the  HR,  the 
HL  layer  had  an  additional  aluminum  reflective  coating  on 
the  a- C  substrate  before  Csl  was  deposited.  Both  types  of 
ACS  had  a  thin  protective  coating,  which  faced  the  optical 
sensor  during  x-ray  exposure.  The  a- C  substrate  side  faced 
the  x-ray  source.  Table  I  is  a  summary  of  the  structures  of  the 
different  Csl  samples. 

B.  PHS  measurements 

PHS  measurements  were  made  on  all  four  Csl  samples. 
The  experimental  apparatus  and  methods  were  the  same  as 
that  described  in  Ref.  19.  Each  Csl  sample  was  placed  di¬ 
rectly  at  the  input  of  the  PMT  (model  R6095,  Hamamatsu 

Table  I.  Summary  of  structural  difference  between  the  four  different  Csl 
samples  used  in  the  PHS  measurement.  Backing  refers  to  the  optical  inter¬ 
face  of  the  Csl  at  the  x-ray  entrance  side.  Front  refers  to  the  side  of  the  Csl 
that  is  in  contact  with  the  optical  sensor,  i.e.,  the  PMT. 


Sample 

Substrate  material 

Backing 

Front  (Sensor  contact) 

FOS/HR 

FOP 

FOP 

FOS/HL 

FOP 

Reflective 

FOP 

ACS/HR 

a-C 

Absorbing 

Transparent  coating 

ACS/HL 

a-C 

Reflective 

Transparent  coating 

Photonics  Corporation),  which  detects  and  amplifies  the  light 
generated  in  Csl  by  each  x-ray  photon.  The  electronic  signal 
at  the  output  of  the  PMT  was  connected  to  a  charge  amplifier 
followed  by  a  pulse-shaping  amplifier.  The  final  shaped  pulse 
had  an  amplitude  proportional  to  the  charge  output  of  the 
PMT.  The  final  output  signal  was  fed  to  a  multichannel  ana¬ 
lyzer  so  that  a  PHS  could  be  obtained.  Monoenergetic  x  rays 
were  used,  which  were  generated  by  K  fluorescence  from 
selected  elements. 

C.  Escape  efficiency  versus  depth 

Here  we  describe  the  method  for  extracting  light  escape 
efficiency  (e)  as  a  function  of  depth  z  from  experimental 
measurements  of  PHS.  With  monoenergetic  x  rays  with  en¬ 
ergy  E  less  than  the  /Ledge  Ek,  the  variations  in  pulse 
heights  observed  in  the  PHS  should  be  due  only  to  the  varia¬ 
tions  in  e  under  the  assumptions  of:  (1)  the  number  of  pho¬ 
tons  m  originally  created  after  absorption  of  an  x  ray  of 
energy  E  is  large  enough  that  the  noise  in  the  photon  creation 
process  can  be  neglected;  and  (2)  the  broadening  of  the  PHS 
due  to  noise  in  the  detection  by  the  PMT  of  the  escaping  and 
collected  photons  is  small  compared  to  broadening  from  the 
variation  in  the  escape  process  itself. 

Assumption  (1)  was  originally  made  by  Swank,7  and  for  a 
screen  in  prompt  emission  mode  and  at  x-ray  energies  such 
that  the  photoelectron  range  is  small  compared  to  the 
thickness,7  should  be  reasonable.  Assumption  (2)  can  be 
tested  by  calculation,  as  in  Drangova  and  Rowlands,  which 
we  outline  here.  As  shown  in  Fig.  1,  the  number  of  electrons 
produced  at  the  photocathode  of  the  PMT  for  a  particular 
pulse  is 

ne  =  mef3r/p,  (5) 

where  e  is  the  escape  efficiency,  /3  is  the  optical  coupling 
efficiency  for  collection  of  emitted  photons  between  the  Csl 
and  the  PMT,  and  rjp  is  the  quantum  efficiency  of  the  PMT 
photocathode.  At  £=30  keV,  taking  the  inherent  energy  re¬ 
quired  per  photon  for  CsI(Tl)  as  —20  eV  (Ref.  19)  results  in 
m  of  —1500.  Assuming  a  mean  escape  efficiency  of  60%, 
and  a  light  collection  efficiency  f3  between  the  Csl  and  the 
PMT  of  80%  results  in  a  number  of  collected  photons  nc  of 
720.  With  a  photocathode  quantum  efficiency  rjp  of  10% 
(Ref.  19)  for  photons  emitted  from  CsI(Tl),  the  number  of 
electrons  produced  would  be  ne=12.  As  will  be  shown  in 
Sec.  IV,  the  mean  number  of  electrons  per  absorbed  x  ray 
finally  detected  by  the  PMT  is  ne~ 70  for  the  ACS  HL 
sample.  The  normalized  standard  deviation  in  the  number  of 
electrons  generated  in  the  photocathode,  due  to  statistical 
effects  in  detection,  is  then  roughly  estimated  as  1  /  \  70 
=  12%.  This  is  small  compared  to  the  spread  observed  in  the 
measured  PHS  for  each  Csl  sample,  which  is  estimated  to  be 
approximately  ±50%.  Although  the  effect  of  PMT  noise  is 
not  completely  negligible,  our  original  assumptions  are  rea¬ 
sonable.  Later  we  will  describe  a  method  for  removing  the 
effect  of  PMT  noise  on  the  PHS  measurements,  so  that  the 
corrected  PHS  can  be  used  to  provide  more  accurate  estimate 
of  the  Swank  factor  and  e. 
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No.  of  electrons  created  by  photo  cathode 

Fig.  4.  Diagram  illustrating  the  integration  of  the  PHS  to  form  the  cumula¬ 
tive  spectrum  C(ne)  for  the  derivation  of  escape  efficiency  vs  depth. 


Given  the  assumptions  above,  the  observed  width  of  PHS 
is  due  to  variations  in  e  versus  depth  z.  If  a  monotonic 
change  in  e  vs  z  is  assumed,  each  value  of  ne  on  the  hori¬ 
zontal  axis  of  the  PHS  corresponds  to  a  unique  value  of 
depth  z.  We  used  the  following  method  to  determine  z(ne), 
i.e.,  z  as  a  function  of  ne.  As  shown  in  Fig.  4,  we  denote  the 
integration  of  the  PHS  for  n  zs  nr  as  C{ne),  which  corre¬ 
sponds  to  the  number  of  x  rays  absorbed  between  the  depths 
of  0 —z  from  the  origin  at  the  exit  surface.  For  the  simplest 
case  of  x  rays  absorbed  uniformly  across  the  thickness  L,  the 
number  of  absorbed  x-rays  C(ne)  is  proportional  to  the  layer 
thickness  z(ne),  and  we  obtain 


z(ne) 


C(ne)L 


(6) 


where  Nq  is  the  total  number  of  x  rays  absorbed  in  the  entire 
phosphor  layer  given  by 

Nq  =  N0{  (7) 


where  N0  is  the  number  of  incident  x-ray  photons  and  /jl  is 
the  energy  dependent  linear  attenuation  coefficient  for  the 
phosphor.  Once  the  depth  z  is  determined  for  a  given  ne 
value,  the  relative  light  escape  efficiency  for  this  depth  can 
be  calculated  using  e=ne/ne{ max),  where  «f(max)  is  the 
maximum  number  of  detected  electrons  at  a  given  x-ray  en¬ 
ergy.  Note  that  the  scale  on  the  e-z  curves  is  not  absolute,  but 
relative  to  the  maximum  signal  obtained,  from  the  side  of  the 
screen  nearest  the  detector. 

We  illustrate  the  above  method  schematically  in  Fig.  5, 
for  clarity.  Suppose  that  the  x-ray  absorption  is  uniform,  and 
that  the  total  number  of  absorption  events  were  N.  Each  ab¬ 
sorption  event  produces  a  large  number  of  photons  which 


travel  along  various  paths  as  illustrated,  but  whose  average 
probability  of  reaching  the  front  surface  and  escaping  is  con¬ 
stant  for  a  given  depth,  and  which  we  define  as  e(z).  Con¬ 
sider  first  the  highest  bin  in  the  measured  PHS,  at  ne= 50 
electrons  and  suppose  it  has  Nl 10  counts.  Then  this  bin  cor¬ 
responds  to  x  rays  that  were  absorbed  within  the  depths  z 
from  0  to  LI  10,  from  the  front  surface.  Taking  ne{ max)  as  the 
number  of  electrons  detected  for  this  highest  bin,  we  plot  e 
=  1,  at  z  of  LI20±L/20.  Consider  next  the  second  highest  bin 
in  the  measured  PHS,  at  ne= 40  electrons  and  suppose  it  has 
3 Nl  10  counts.  Then  this  bin  corresponds  to  x  rays  that  were 
absorbed  between  the  depths  of  LI  10  and  4Z./10  from  the 
front  surface.  We  plot  r'  =  0.8,  at  z  of  5 L /  20  ±  3  L/  20, and  simi¬ 
larly  for  the  rest  of  the  bins  of  the  PHS  histogram. 

We  can  extend  the  above  analysis  to  nonuniform  (expo¬ 
nential)  x-ray  absorption.  For  the  front  screen  geometry 
shown  in  Fig.  1,  (so  called  because  the  photons  are  emitted 
from  the  side  opposite  to  where  the  x  rays  enter),  the  C{ne)  is 
given  by 

C(ne)=N0[e-^-e^L].  (8) 


Substituting  Eq.  (7)  into  Eq.  (8)  we  obtain 

CM‘  a • 


(9) 


The  value  of  Nq  can  be  determined  by  integrating  the  entire 
PHS,  which  gives  the  total  number  of  x  rays  detected  in  the 
experiment.  Solving  Eq.  (9)  for  z,  we  obtain 


z(ne)  =  —  In 
U- 


1  + 


C(ne)(^L~  1) 
Nn 


(10) 


The  following  steps  are  used  to  derive  c  as  a  function  of  z 
from  the  measured  PHS:  (1)  determine  /jl  for  the  x-ray  en¬ 
ergy  used  in  the  PHS  measurement;  (2)  calculate  Nq  by  in¬ 
tegrating  the  entire  PHS  to  obtain  the  total  number  of  x  rays 
absorbed;  (3)  choose  a  value  for  ne  on  the  horizontal  axis  of 
PHS,  and  calculate  C{ne)  by  integrating  the  PHS  for 
and  (4)  calculate  depth  z  for  the  given  ne,  using  Eq.  (10)  and 
determine  e  for  this  depth  using  e  =  ne/ne{ max). 

The  above  analysis  assumed  that  the  effect  of  PMT  detec¬ 
tor  noise  on  the  PHS  measurements  was  negligible.  Now  we 
describe  a  method  for  estimating  and  removing  the  effect  of 
PMT  detector  noise  from  the  measured  PHS,  so  that  the 
e-z  curves  can  be  extracted  more  accurately.  As  shown  in 
Fig.  1,  the  mean  number  of  photons  collected  by  the  PMT  is 
nc.  Assuming  that  the  input  photon  stream  is  noiseless,  the 


EE 


Relative  optical  efficiency 


Fig.  5.  Schematic  illustrating  the  construction  of  the 
relative  optical  efficiency  via  the  measured  PHS. 
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number  of  photoelectrons  created  in  the  PMT,  ne,  will  vary 
due  to  the  binomial  distribution  in  the  quantum  absorption 
process  of  the  photocathode,  with  efficiency  r/p.  The  variance 
in  the  number  of  photoelectrons  crne  is  then  given  by" 

vie  =  VpO  -  Vp)nc-  (11) 

Further  noise  is  added  due  to  the  statistical  fluctuation  in  the 
gain  of  the  PMT  dynode  chain,  which  is  usually  quantified  as 
a  multiplicative  excess  noise  factor  (NF)  which  multiplies 
the  total  variance  given  above.  NF  is  dominated  by  the  first 

dynode  stage  of  the  PMT,  and  can  be  estimated,  for  Poisson 

20 

dynode  statistics,  from  the  first  dynode  gain  S  using 

NF  ££  S/(S-  1) .  (12) 

For  the  PMT  and  power  supply  voltage  used  in  our  experi¬ 
ments,  S  was  estimated  to  be  3.76,  resulting  in  a  NF  value  of 
—  1.36."1  The  measured  PHS,  and  the  estimated  PMT  noise 
were  both  fitted  to  Gaussian  functions.  A  corrected  PHS  was 
derived  by  deconvolution  of  the  PMT  noise.  The  e(z)  extrac¬ 
tion  procedure,  as  described  above,  was  then  repeated  using 
the  corrected  PHS. 


D.  Calculating  /0PD  from  escape  efficiency 

The  most  straightforward  method  for  obtaining  the  optical 
factor  /OPD  is  from  the  moments  of  the  measured  PHS  using 
Eq.  (3)  for  E<EK,  where  IAKD  can  be  assumed  to  be  unity. 
Also,  for  energies  E>EK,  the  factor  7AED  can  be  calculated 
by  simulation  of  characteristic  x-ray  emission  and  reabsorp- 
tion  processes,  and  /0PD  can  then  be  estimated  indirectly 
using  Eq.  (4).  However  in  screen  optimization  procedure,  we 
sometimes  only  have  the  estimated  light  escape  efficiency 
versus  depth  based  on  a  particular  screen  design  model.  We 
will  show  here  that  /0PD  can  also  be  obtained  from  the  e-z 
curve.  This  will  provide  a  theoretical  means  to  determine  the 
Swank  factor  for  different  screen  designs  without  the  need 
for  explicitly  calculating  the  PHS. 

The  amplitude  of  the  scintillation  pulse,  i.e.,  the  number 
of  escaping  light  photons  for  each  absorbed  x  ray,  n(z)  is 
given  by: 

n{z)  =  m(E)e(z),  (13) 

where  m{E )  is  the  number  of  luminescent  photons  produced 
by  absorption  of  an  x  ray  of  energy  E  in  the  material,  and 
e(z)  is  the  escape  efficiency  from  depth  z.  The  /th  moment  of 
the  optical  scintillation  spectrum  is  then 

Mj  =  J  Nx(z)mJe(zydz,  (14) 

where  Nx(z)  is  the  number  of  x  rays  absorbed  at  depth  z. 
Neglecting  the  nonuniformity  in  Nx{z),  we  obtain 

M0  =  Nxf  dz,  (15) 
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M]  =  Nxm 
and 

M2  =  Nxm 


Then  the  optical  factor  becomes 

M\  1  /L2Je(z)dz2 
OPD“m2M0~  1  /Lje(z)2dz  ' 


(16) 


(17) 


(18) 


Equation  (18)  indicates  that  70Pd  (in  the  case  of  uniform 
x-ray  absorption  as  a  function  of  depth  z)  is  determined  by 
the  square  of  the  depth-average  escape  efficiency  divided  by 
the  average  square  escape  efficiency 

IOPD  =  (e)2/(e2).  (19) 


To  demonstrate  the  derivation  of  70Pd  from  e(z),  we  used 
several  simple  examples  of  screen  optics  design:  (1)  White 
backing  turbid  screen:  100%  reflectance  for  the  screen  back¬ 
ing,  no  bulk  absorption  and  no  internal  reflection  at  exit;  (2) 
Black  backing  turbid  screen:  100%  absorption  for  the  screen 
backing,  no  bulk  absorption  and  no  internal  reflection  at  exit; 
(3)  Linear  escape  efficiency:  e(z)  decreasing  linearly  with 
depth,  with  a  mean  value  of  e  and  total  variation  2Ae,  where 
e(z)  is  given  by 

e(z)  =  e-  Ae|— -lj  (20) 


and  (4),  Turbid  screen  with  dye  absorption  added:  the  case  of 
added  bulk  absorption  illustrated  in  Fig.  2(a). 


IV.  RESULTS  AND  DISCUSSION 
A.  PHS  measurements 

PHS  measurements  were  performed  for  a  range  of  inci¬ 
dent  x-ray  energies  between  25  and  50  keV.  Shown  in  Fig.  6 
are  representative  PHS  measurements  for  the  ACS  HL 
sample  at  x-ray  energies  of  25,  30,  40,  43,  46,  and  49  keV. 
Figure  6  shows  that  at  E>EK  of  Cs  and  I  (36  and  33  keV, 
respectively),  there  is  a  second  peak  in  the  PHS  due  to  the 
escape  of  TGfluorescence  photons.  Since  the  purpose  of  this 
article  is  to  investigate  7OPd  and  extract  the  light  escape  ef¬ 
ficiency  for  different  screen  designs,  only  the  PHS  for  E 
<Ek  were  used.  The  effect  of  K  fluorescence  is  avoided 
which  simplifies  the  data  analysis.  The  measured  PHS  for  all 
four  Csl  samples  are  shown  in  Fig.  7  for  E=25  and  31  keV. 
The  PHS  are  plotted  as  the  number  of  counts  P{ne)  versus 
the  size  of  the  scintillation  pulse  ne  (given  in  number  of 
electrons  detected  by  the  PMT).  The  FOS  samples  were 
coupled  to  the  PMT  through  the  FOP  substrate,  which  had  an 
optical  transmittance  of  60%.  The  ACS  samples  were 
coupled  to  the  PMT  through  a  transparent  protective  coating, 
which  had  minimal  attenuation  of  light.  To  allow  for  easier 
comparison  between  the  FOS  and  ACS  sample  types,  the 
PHS  for  the  FOS  samples  (HL  and  HR)  are  replotted  in  Fig. 
8  after  correction  for  the  60%  FOP  transmittance.  Figure  8 
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Fig.  6.  Measured  PHS  for  the  ACS  HL  sample  at  six 
different  x-ray  energies. 


shows  that,  after  correction,  the  ACS  HL  sample  produced 
the  most  light  photons  per  x  ray,  and  FOS  HR  and  ACS  HR 
the  least.  Although  the  mean  light  output  from  the  ACS  HR 
is  comparable  to  that  from  the  FOS  HR,  it  has  a  broader 
distribution  in  the  PHS.  Given  that  the  PHS  measurements 
for  all  four  samples  were  made  using  the  same  apparatus, 
and  the  samples  were  made  with  the  same  scintillating  ma¬ 
terial,  it  can  be  concluded  even  before  detailed  analysis  that 
the  magnitude  of  e  was  markedly  different  between  screen 
types. 

B.  Escape  efficiency  versus  depth 
1.  Without  PMT  noise  correction 

The  PHS  measurements  of  the  four  Csl  samples  were 
used  to  obtain  e-z  curves.  The  results  are  shown  in  Fig.  9,  for 
E= 25  and  31  keV.  Figure  9  shows  that  the  e-z  curves  de¬ 
rived  from  the  two  sets  of  PHS  measurements  at  E= 25  and 
31  keV  are  nearly  identical,  supporting  the  theory  that  the 
optical  escape  process  is  independent  of  the  incident  x-ray 
energy.  The  value  of  e  decreases  approximately  linearly  with 
Z  for  all  the  Csl  samples.  There  is  some  bending  in  the  curves 
at  the  high  and  low  ends  of  z,  which  correspond  to  the  tails 
of  the  PHS.  In  an  ideal  system  with  no  photon  noise  and  no 
measurement  noise,  an  e(z)  which  was  perfectly  linear  with 
depth  would  result  in  a  PHS  with  a  flat  top,  and  vertical 


Fig.  7.  (a)  Measured  PHS  for  four  different  Csl  samples  at  £=25  keV.  (b) 
Measured  PHS  for  four  different  Csl  samples  at  £'=31  keV. 
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(b)  Number  of  Electrons  ne 

Fig.  8.  (a)  Measured  PHS  for  four  different  Csl  samples,  with  results  for  the 
FOS  samples  adjusted  to  account  for  optical  transmittance  of  FOP  substrate 
at  £= 25  keV.  (b)  Measured  PHS  for  four  different  Csl  samples,  with  results 
for  the  FOS  samples  adjusted  to  account  for  optical  transmittance  of  FOP 
substrate  at  £=31  keV. 

edges.  In  reality  noise  sources  exist,  resulting  in  unsharpness 
in  the  PHS  edges  and  bending  in  the  derived  e-z  curves.  In 
the  next  section  the  effect  of  measurement  noise  on  the  ob¬ 
served  PHS  and  derived  e-z  curves  is  addressed. 

It  should  also  be  understood  that  the  relative  optical  effi¬ 
ciency  plotted  in  Fig.  9  is  not  an  absolute  number,  but  a 
fraction  relative  to  the  highest  experimentally  observed  num¬ 
ber  of  detected  electrons,  ne(max).  There  is  some  uncertainty 
(10%)  in  picking  nf(max)  because  of  the  unsharpness  in  the 
edges  of  the  measured  PHS. 

2.  With  PMT  noise  correction 

It  should  be  noted  that  the  uncorrected  PHS  measure¬ 
ments  shown  in  Fig.  7  were  used  to  estimate  the  effect  of 
PMT  detector  noise.  Figure  10  shows  the  comparison  be¬ 
tween  the  measured  PHS  of  the  ACS  HL  sample  and  its 
Gaussian  fitting,  which  is  in  reasonable  agreement  with  the 
measurement  and  facilitates  the  deconvolution  process  for 
removing  the  effect  of  PMT  noise  from  the  PHS  measure¬ 
ments.  Also  shown  in  Fig.  10  is  the  PMT  noise  corrected 
PHS,  which  is  narrower  than  the  original  PHS  measurement 
data.  The  e(z)  derived  from  the  corrected  PHS  are  shown  in 
Fig.  9(c)  for  E=  25  keV.  It  can  be  seen  that  the  bending  in 
the  curves  at  the  high  and  low  ends  of  z,  as  well  as  the  slope 
of  the  curves  in  the  linear  region,  are  both  reduced  compared 
to  the  uncorrected  results  in  Fig.  9(a). 

Screen  optics  models  provide  a  link  between  e(z)  and 
optical  design  parameters,  e.g.,  substrate  reflectance  or  bulk 

o 

absorption.  Such  models  have  been  useful  for  the  optimiza- 
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tion  of  systems  employing  powdered  phosphor  screens. 
Rowlands  and  Yorkston"  have  shown  in  arguments  supported 


Fig.  9.  Light  escape  efficiency  vs  depth  derived  from  the  measured  PHS  for 
four  Csl  samples  at:  (a)  £=25  keV;  (b)  £=31  keV;  and  (c)  £=25  keV,  with 
PHS  corrected  for  the  effect  of  PMT  noise. 


by  experimental  data  that,  although  Csl  screens  exhibit  better 
spatial  resolution  than  powder  screens  of  the  same  thickness, 
they  show  an  optical  behavior  intermediate  between  a  perfect 
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Fig.  10.  Comparison  of  PHS  with  and  without  PMT  noise  deconvolution. 
The  data  points  are  the  measured  PHS  for  the  ACS  HL  sample.  The  solid 
line  is  a  Gaussian  fit  to  the  PHS  data,  and  the  dashed  line  (with  smaller 
width)  is  the  PHS  corrected  for  PMT  noise. 


Fig.  11.  Schematic  illustration  of  the  escape  efficiency  vs  depth  for  three 
simple  cases:  (a)  100%  backing  reflectance,  no  bulk  absorption  or  internal 
reflection  at  exit;  (b)  100%  backing  absorption,  no  bulk  absorption  or  inter¬ 
nal  reflection  at  exit;  and  (c)  escape  efficiency  decreases  linearly  vs  depth. 


fiber  optic  and  a  turbid  screen,  with  MTF  decreasing  with 
increasing  thickness.  Although  it  would  be  helpful  to  develop 
analytic  physical  models  for  anisotropic  (columnar  struc¬ 
tured)  screens  like  Csl,  we  may,  with  caution,  use  insights 
from  the  isotropic  powder  models  to  help  interpret  the  de¬ 
rived  e-z  curves.  For  example,  the  e(z)  for  ACS  HR  in  Fig.  9 
above  drops  to  near  zero  at  the  backing,  whereas  ACS  HL 
has  a  generally  higher  e(z)  in  the  bulk,  approaching  a  non¬ 
zero  value,  for  z  near  the  back.  These  are  the  kinds  of  be¬ 
havior  expected  (as  in  Fig.  2)  for  screens  comprising  phos¬ 
phors  on  nonreflecting  and  reflecting  substrates,  respectively. 
The  model  interpretation  is  in  agreement  with  the  optical 
design  of  the  ACS  HR  and  ACS  HL  samples,  which  are 
optimized  for  high  resolution  and  high  light  output,  respec¬ 
tively.  If  we  compare  the  two  screens  optimized  for  high 
resolution  we  find  that  compared  to  the  ACS  HR,  and  also  in 
contrast  to  the  predicted  e(z )  behavior  for  a  turbid  phosphor 
HR  screen  with  a  nonreflective  backing,  the  FOS  HR  has  an 
e(z)  that  varies  much  less  versus  z.  The  substantially  lower 
e(z )  of  the  ACS  HR  near  the  back  is  because  the  Csl  layer 
was  deposited  on  a  black  substrate  (a-C),  which  has  much 
higher  optical  absorption.  The  reduced  z  dependence  of  e(z) 
for  the  structured  phosphor  FOS  HR  screen,  compared  to  the 
strong  depth  dependence  expected  for  powder  screens  with 
nonreflective  backings  shown  in  Fig.  2(b),  is  in  agreement 
with  predictions  based  on  Monte  Carlo  calculations  for  light 
scattering  in  homogeneous  compared  to  columnar  screen 
structures.  This  behavior  is  due  to  partial  waveguiding  of 
light  up  and  down  the  columns. 

C.  Calculating  /0PD  from  escape  efficiency 

The  e-z  curves  for  the  first  three  simple  cases  listed  in 
Sec.  HID:  (1)  100%  reflectance  for  the  screen  backing,  no 
bulk  absorption,  no  internal  reflection  at  exit;  (2)  100%  ab¬ 
sorption  for  the  screen  backing,  no  bulk  absorption,  no  inter¬ 
nal  reflection  at  exit;  and  (3)  e(z)  decreasing  linearly  with  z. 


are  shown  in  Fig.  11.  The  e(z)  for  the  first  case  is  unity 
across  the  full  range  of  z,  resulting  in  /0pd=  I .  The  e(z)  for 
the  second  case  decreases  from  unity  at  z= 0  to  zero  at  the 
back  of  the  screen,  i.e.,  z=L,  and  we  obtain  7OPD=0.75,  as 

g 

originally  pointed  out  by  Swank.  For  the  third  case,  /0PD  can 
be  derived  using  Eqs.  (19)  and  (20)  as 

1/^opd  =  1  +  ( l/3)(Ae/e)2.  (21) 

For  the  fourth  simple  case  of  added  bulk  absorption  illus¬ 
trated  in  Fig.  2(a),  the  calculated  /0PD  decreases  from  0.75  to 
0.30  with  increasing  bulk  absorption  for  the  examples 
shown.  By  adding  even  more  bulk  absorption,  the  Swank 
factor  can  eventually  be  reduced  to  zero  in  the  limit  of  high 
optical  absorption. 

For  the  linear  case,  as  shown  in  Eq.  (21)  above,  70Pd 
depends  only  on  the  fractional  variation  in  e(z).  Thus  in  this 
case,  either  improving  e  at  a  constant  Ae,  or  reducing  Ae  at 
a  constant  e,  will  both  improve  70PD.  The  experimental  Csl 
screens  above  have  derived  e-z  curves  that  resemble  the  lin¬ 
ear  case.  For  the  results  in  Fig.  9,  the  e-z  curves  were  fitted 
with  a  straight  line,  as  shown  for  example  in  Fig.  9(b),  and 
the  70pd  from  Eq.  (21)  was  compared  with  that  obtained  in 
the  usual  way  from  the  moments  of  the  PHS.  The  results 
were  generally  equivalent,  as  shown  in  Table  II. 

The  70Pd  was  also  calculated  for  each  sample,  using  the 
PHS  that  were  corrected  for  the  PMT  detector  noise  (as  de¬ 
scribed  in  Secs.  Ill  C  and  IV  B).  The  results  are  shown  in 
Table  II.  Neglecting  the  detector  noise  underestimates  70PD 
by  5%-10%  for  the  FOS  samples,  and  by  <3%  for  the  ACS 
samples,  which  had  a  wider  PHS  distribution  than  the  FOS 
samples.  From  the  e-z  curves  in  Fig.  9(c),  it  can  be  seen  that 
the  FOS  HR  and  FOS  HL  have  essentially  the  same  Ae,  and 
the  higher  e  for  the  FOS  HL  sample  leads  to  a  slightly  higher 
70Pd-  Although  the  ACS  HL  has  a  slightly  larger  Ae  com¬ 
pared  to  the  FOS  samples,  it  has  an  70Pd  value  comparable  to 
that  of  the  FOS  samples  because  it  has  the  highest  e  for  all 
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Table  II.  Comparison  of  optical  excess  noise  factors  obtained  from  measured  Pulse  Height  Spectra  for  CsI:T  1 
samples  of  various  optical  designs,  with  the  noise  factors  obtained  from  (i)  measured  uncorrected  Pulse  Height 
Spectra,  (ii)  linear  tit  to  uncorrected  derived  e-z  curves,  (iii)  Pulse  Height  Spectra  corrected  for  the  effect  of 
PMT  detector  noise. 


Sample 

From  uncorrected  PHS 

From  linear  fit  to 
(uncorrected)  e-z  curve 

From  corrected  PHS 

FOS/HR 

0.90 

0.91 

0.97 

FOS/HL 

0.93 

0.94 

0.98 

ACS/HR 

0.86 

0.82 

0.86 

ACS/HL 

0.94 

0.95 

0.97 

the  samples  used  in  our  experiments.  The  ACS  HR  sample, 
which  has  the  highest  Ae  due  to  its  absorptive  backing,  has 
the  lowest  /OPD. 

The  goal  of  our  investigation  is  to  understand  the  inter¬ 
play  of  screen  parameters  so  as  to  permit  the  optimization  of 
the  performance  of  Csl  layers  for  different  DR  detector  ap¬ 
plications.  It  can  be  seen  from  Eq.  (2)  that  to  optimize  the 
DQE(O),  we  need  to  maximize  both  the  quantum  absorption 
efficiency  77  and  the  Swank  factor  As.  Equation  (20)  shows 
that  to  maximize  I0PD ,  we  need  to  maximize  the  mean  escape 
efficiency  e,  and  minimize  its  variation  Ae.  Screen  optics 
models  indicate  that  one  way  to  accomplish  this  is  to  use  a 
screen  with  high  quantum  gain  m,  minimal  bulk  absorption 
and  a  highly  reflective  backing.  For  the  four  samples  used  in 
our  experiment,  the  two  ACS  samples  have  larger  Ae,  corre¬ 
sponding  to  a  larger  variation  of  e  with  depth.  The  ACS  HL, 
FOS  HL,  and  FOS  HR  samples  all  provide  very  high  levels 
of  70pd  and  A s. 

It  is  noted  that  the  /OPD  values  obtained  for  the  Csl 
samples  are  quite  high,  i.e.,  3=0.9  except  for  ACS  HR.  These 
values  compare  favorably  with  the  range  of  I0PD  of  0.6  -  0.9 
for  good  quality  powder  phosphor  screens.18  The  higher  /0PD 
for  Csl  is  most  likely  connected  with  the  light  channeling 
properties  of  the  Csl  needle  structure,  which  evidently  re¬ 
duces  the  variation  of  escape  efficiency  with  depth.  This  may 
be  due  to  reduced  optical  path  length  and  bulk  absorption  for 
photons  produced  at  some  depth  from  the  exit  surface  be¬ 
cause  their  travel  is  at  least  partially  restricted  to  columnar 
paths. 

To  optimize  the  performance  of  an  imaging  detector,  in 
which  performance  at  finite  spatial  frequencies  is  important, 
we  must  go  further  than  DQE(0)  and  consider  the  DQE  as  a 
function  of  spatial  frequency.  This  will  be  the  subject  of 
future  work. 

V.  CONCLUSIONS 

A  theoretical  model  was  developed  for  the  light  escape 
efficiency  of  x-ray  scintillators,  and  a  method  for  deriving 
e(z )  from  PHS  was  described.  The  method  was  applied  to 
experimental  PHS  measurements  from  Csl  samples  with  four 
different  designs.  The  results  showed  that  e(z)  varies  essen¬ 
tially  linearly  as  a  function  of  z  in  the  Csl  samples,  and  that 
the  magnitude  of  variation  is  relatively  small,  leading  to  very 
high  values  for  both  the  /0PD  and  Swank  factor.  The  effect  of 


PMT  noise  on  PHS  measurement  was  estimated  and  found  to 
result  in  a  small  but  noticeable  underestimation  (—5%)  of 
the  Swank  factor. 
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